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Abstract It is shown how an explicit solution to the problem of finding the minimum holding 
time and best temperature profile may be obtained for a tubular or batch reactor. The form 
which this takes allows a complete parametric study for given kinetics and full results are pre- 
sented for a first-order reversible reaction. For the sequence of stirred tank reactors both 
the minimum holding time and optimum temperature of each stage can be found by a 
simple graphical construction. 


Résumé—-Le temps de séjour minimum et le meilleur profil de températures peuvent étre obtenus 
par une solution mathématique explicite pour un réacteur tubulaire ou discontinu. Cette solution 
permet une étude parameétrique compléte pour des cinétiques données, et des résultats complets 
sont donnés pour une réaction reversible du 1 er ordre. Pour la série de réacteurs & cuves agités 
une construction graphique simple permet de trouver le temps de séjour nimimum et la température 
optimum pour chaque stade. 


Zusammenfassung —Es wird gezeigt, wie man eine explizite Losung fiir die Ermittlung der 
kleinsten Verweilzeit und des giinstigsten Temperaturprofils fiir einen Rohrreaktor oder einen 
diskontinuierlichen Reaktor erhalt. Hieraus kénnen bei gegebener Kinetik alle Parameter 
untersucht werden. Es werden die Ergebnisse fiir eine reversible Reaktion I Ordnung angegeben. 
Fiir eine Riihrkesselkaskade kann die kiirzeste Verweilzeit und die optimale Temperatur fiir 
jeden Kessel durch eine einfache graphische Methode ermittelt werden. 


Ir micut be felt that little more could be added to give the optimum choice both of temperature 
to the discussion of optimum temperature and holding time at each stage. The example of a 
profiles in tubular or stirred tank reactors, first-order reversible reaction has been used to 


when only a single reaction is taking place. 
It has long been known that the temperature 
should always be chosen to make the reaction 
rate as large as possible at each point or stage 
[1-8], and the consequent calculations are not 
difficult. Two observations, however, prompt 
a further remark on the subject. The first is 
that an explicit solution in terms of easily calcul- 
able integrals can be found and this allows a full 
parametric study of any particular reaction ; 
the second, that Denpicn’s notion of the 
maximization of rectangles [4] can be adapted 


illustrate the methods, but these are perfectly 
general and the explicit solution applies to the 
kinetic expression commonly given for a homo- 
geneous reaction. 


1. Tue Reaction anp Reactor Systems 


We consider a single reaction between n 
chemical species A,, A,, ..., 4,, which we write 
in the form 0 = x, A; and observe the conven- 


tion that the stoichiometric coefficients, x, of 


products should be positive whilst those of 
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reactants are negative. If ¢; denotes the concen- 
tration (moles unit volume) of 4,;. and ey its 
concentration in the feed (or in some reference 
mixture from which the feed might be derived by 
reaction to a certain extent), then 


= Ce + (1) 


where ¢ is a measure of the extent of the reaction. 
The rate of the reaction, which we assume takes 
place without change of volume, is by definition 
de dt and is a function of all the concentrations 
and of the temperature. We therefore write 


rie. T). (2) 


dt 


The expression commonly used to fit homogeneous 
reaction rate data is 


r(e, T) (T) ky(T) ii (3) 


where 8, and y, are experimentally determined 
indices with more or less theoretical backing, 
but in general restricted only by the equilibrium 
consistency condition, y, — 8B, =~ x. The explicit 
solution we give will be restricted to the kinetic 
expression (3). but the same method is applicable 
to any form of r(c, T). 

Equation (2) describes the course of the 
reaction for a batch or tubular reactor; for the 
batch reactor ¢ is the actual time, for the tubular 
it is the length of reactor up to a given point 
divided by the linear velocity of flow. The 
optimum problem is to choose T as a function 
of ¢ so that cither the time (or reactor length) 
required for a given conversion is least, or the 
conversion is a given time is greatest, AMUNDsoN 
and Bitous have shown the equivalence of these 
problems [5] and we shall consider the first 
formulation. The admissible choice of tempera- 
ture will usually be restricted by upper and lower 


bounds, 
T, s T > T* (4) 


The sequence of stirred tank reactors is 
described by a set of difference equations analo- 
gous to (2). It is convenient to number the stages 
backwards, the first reactor of a battery of R 
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reactors being labelled reactor R and the last, 
reactor 1. The state of the feed to reactor r 
is then specified by ¢,,,. the extent of reaction in 
the previous reactor, whilst the composition 
within it, and so also its product composition, is 


denoted by ¢,. (Fig. 1 illustrates this). The 


| 
R 
—e-- 


Cc 


Con | 


Fic. 1. The stirred tank sequence. 

design of the sequence involves the choice of 
the holding times @,....@, and operating tem- 
peratures 7, ... so that a given conversion 
may be attained with minimum total holding 


time O, Se, A mass balance on the 
reactor gives the difference equation 
Coty Ce — T,) (5) 
or 


In the latter form we see that, if @, is allowed to 
tend to zero, the differential equation (2) arises 
as the limit of the difference equation (5). Solv- 
ing (5) for 0, we have 


= as 1 (6) 


and given the states of the feed and product, 
Cpe, and ¢,, we have to choose the temperatures 
T,,... Ty and the extents of reaction cy. ... 
which will make 9, a minimum. This involves 
the choice of 2R — 1 quantities, which concurs 
with the other formulation of the problem 
(namely the maximization (¢, — ¢g.,) for given 
O,). for in this case the R temperatures have to 
be chosen and R — 1 of the holding times, the 
last being fixed by the value of @,. 


For completeness we shall give a brief demon- 
stration of the fact that 7 should always be 
chosen to maximize r(c, 7). Earlier published 
proofs have not taken into account the restriction 
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(4) and, though true, it is not obvious that this 
does not change the result. For the batch reactor 
we observe that we have to minimize. 


J T) 

by choice of T(c) in the range ¢ 
(Here we denote the product state by ¢, as in 
the case of the stirred tank sequence, but the 
feed state by ey). T as a function of ¢ will then 
be obtained by integrating equation (2) to get 
the relation between ¢ and t. ¢, being fixed the 
minimum value of this integral will depend only 
on ¢, and we may write 


fe) — Min (7) 

Now the range of integration may be broken 
up into a short interval (¢, y) and the remainder 
(y. ¢). Moreover it is clear that whatever choice 
of T (ce) we make in the first interval we cannot 
do better than to make the best choice in the 
remainder of the interval to which the feed 
condition is ¢ = y. Indeed by trying different 
values of T (c) in (ey, y) and always using the 
optimum choice in (y, ¢,) we shall eventually 
hit on the optimum choice for the first interval 

also. Thus 


de 


Now let y» > cy and let 7; be the inlet tempera- 
ture. Then 


Mite 
fey) = Min | + 


+(y —&) Sf (gq) + Oly — 


-f (ey) Min r(c. T)) 
Since <0 this means that 7, must be 
chosen to make r(cy, Ty) as large as possible. 
This demonstration is valid under the restrictions 
(4), whereas proofs involving differentiation break 
down under this restriction. 


For the stirred tank sequence we again argue 
that when the optimum choice has been made the 
minimum holding time will be a function of the 
feed state only. Let. 


Sr (Cris) = Min Op (8) 


then whatever choice of cg, and T'y may be made 
in the first reactor we cannot do better than use 
the optimum policy for the remaining (r — 1) 
reactors, with respect to their feed cp. Thus 


Sn Min + Su (en)| (9) 


Now cz and Tz are independent choices and the 
second term in the bracket depends only on 
Cy It is thus clear that so far as the choice of 
Tp is concerned it should always be made to 
maximize r(c, 7). This may be determined by 
any means or be known experimentally ; how- 
ever, with the kinetic expression (3) we may 
find an explicit solution. 


3. Ture Maximum Rate ror Given 
COMPOSITION 


In the general reaction-rate expression (3) 
k, and k, are assumed to have the Arrhenius 
form kj = k* exp (— E,/RT), i= 1,2. Then, 
since dk,/dT — E,k,/ RT*, the maximum value 
of r is for T satisfying 

F, ky ~ Eyk, i = 0 (10) 
provided this falls in the interval 7',, 7* and 
makes the second derivative of r negative. But 
the second derivative of r at this extremum is 
proportional to 


it — k, it 


For an exothermic reaction when r, given by (3), is 
positive and the expression in (10) is zero this is 
certainly negative. For an endothermic reaction 
with r positive, r always increases with increase 
of T so that T should always have its maximum 
value T*. The case of the endothermic reaction 
is thus rather trivial. 

Equation (10) may be solved for T in the form 
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ky _ 

since y, — Bj ~ x. Let us denote the dimen- 

sionless temperature by + = RT /(— AH), then 

since AH = E, — E,. 

— _RT = 
(— AM) 


exp (E, — E,) RT — 


[in (11) 


The maximum reaction rate may now be written 


Ps, 
R (ce) = ke (12) 
where 


(13) 


The concentrations ¢; are of course given by (1) 
and it is convenient to make ¢ dimensionless by 
the following substitutions. Let 


1 


then 


C/Ce 


E,) 
R (0) = ATT (uy + (12) 
where 


(p + k,*? 


py-l (13) 


has the dimensions of (time)'. The value of + 
can now be written as 


(14) 


‘= 


where 
It thus appears that for reaction mixtures of 
given proportions (,,) and for given kinetics 
(8,;, y) there are only two parameters p and « 
that need to be varied to study the whole class 
of optimum problems, for A may be used to 
render time dimensionless, 
If the value of T given by (11) does not lie in 
the interval (7',, T*) then T must take one of 


(15) 


these values. In this case R(c) the maximum 
value of r is either R* or R, given by (3) with 
T = T* or T, respectively. 

For the case of a first-order reversible reaction 
O0=4,—Ay 1 B= =O, the 
stoichiometric relations imply ¢, + + 
In this case it convenient to put 
+ 20) = (Cre + + and we 
have 

R(c) = 


= k*, Ey/k,* BE, J 


4. Tue Temperature Prorite 
From the relation (14) it follows that 


dr 

is always negative. Hence there will be two 
definite values 7, and r* at which + — +* and 
*(T — T*, T,) respectively, and only within 
the interval (2,, .7*) will the expression (12) be 
used, Thus using the optimum policy and letting 
At = 2, equation (2) becomes 


x," 


= (2), 2 < (17) 


(18) 


where = R(c)/eg A. 


From this the solution is immediately obtained 
by quadratures, and z, the dimensionless time or 
reactor length, is the area under the curve 1/P 
shown in Fig. 5 between the ordinates with 
values of 2 corresponding to ¢, and ¢,. The 
temperature profile may be immediately found 
by graduating this curve in values of + calculated 
by (14). It is however interesting to see how a 
complete parametric study of a given reaction 
may be made and for this we may consider the 
first-order reversible reaction. 


13 
dx 
— 
dr 
P(z), % <2 <2", 
dr 
= Pele), <x, (19) 
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Fic. 2. Minimum length of reactor for first-order reactions. 


Disregarding for the moment the restrictions on 
temperature we have from equations (16) and 
(18) 


[va 


the limits of the integration being the values of 
= €,/(Cy + Cy) in the feed and _ product. 
This integral depends only on p and hence a 


whole family of curves can be drawn as in Fig. 2 
giving the relation between z and the upper 
limit of integration. Since p = E,/(E, — E,) is 
positive for an exothermic reaction, the curves 
are as shown, sloping upward from the origin 
y — 0 to the asymptote y = 1. 

In Fig. 3 the temperature policy 7 is shown 
for the same range of y. Here the only parameter 
is x and we have a family of curves all tending to 


o4 


os 


mole fraction, y 
Fic. 3. Optimum temperature profiles for first-order reactions. 


Chem, Engng. Sci. Vol. 18, No. 4. June, 1961. 


201 


é 
— 
13 mole fraction, y 
26 
fj 
8 = + 
| 
| | 
3 
2 I | 
02 03 = = OF 08 ar) 
|_| 


R. Aris 


mole fraction, y 


Fic. 4. Construction of the optimum temperature 
gradient. 


7=0 as y=1 and having asymptotes at 
y—(1 + Curves from the two figures are 
combined in Fig. 4. From this it is immediately 
clear what is the range of validity of (18). Given 
7* we find the point A on the appropriate « 
curve and this gives y,; likewise y* is obtained 


ke 


Fic. 5. Construction for the holding times of a 
sequence of stirred tanks. 


from rt, via the point C. Thus the relevant part 
of the z curve for the appropriate p is the arc 
BD. Outside this range we may immediately 
integrate either equation (17) or (19) to give 


A jk, — (hy + 
where k, and k, are evaluated either at tempera- 
ture T = T, or at T — T*. 

The same treatment can be accorded the more 
general reaction, where the solution of (18) now is 


[ius + dar, 


an integral presenting little computational 
difficulty. We notice that for exothermic reac- 
tions 6; px; — 8; will generally be negative for 
reactants and positive for products. The general 
form of the curve is as before with = tending to 
infinity as 2 approaches the least value of 
— p/x, Of course this point is never reached 
for the optimum temperature here would be 
absolute zero, so that the lower bound 7, inter- 
venes and beyond this the best that can be done 
is to allow the reactants to approach equilibrium, 


5. Tur Srirrep Tank Sequence 


The functional equation (9) by which the 
minimum holding time fz (¢g, ,) may be calculated 
has already been derived and since we now know 
that r should have its maximum value R we may 
write 

. — Ces 
(Cr+) Min ats 

R (cp) 
where the minimization is by choice of cg in the 
range Cas, — To make this dimension- 


less let A O, and dp Min dx. Then 
1 


(en) 


[tp — Tarn 
$n Min P 
The holding time ¢, of each stage is thus the 
area of the rectangle standing on the base 
and of height 1/P(2,). For any choice 
of 2, @-4.+.+@y these rectangles can be drawn 
as in Fig. 5 and serve to evaluate the holding 


+ (20) 
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time of each stage. It is clear that the total 
holding time for the sequence of stirred tanks 
must exceed that of the tubular reactor, which is 
the area under the curve itself, and the problem 
is to minimize this area by the correct choice of 
This idea of minimization of rectangles 
derives from Dennicu’s recent presentation of a 
method of maximization of rectangles [4], which 
applies to the equivalent problem of maximizing 
the yield for given holding time. Denpicn is 
concerned with more complicated reaction systems 
under isothermal conditions. For R — 1 there 
is no choice, for to attain a given conversion a 
holding time @, (a, — ®)/P(a,) is required. 
We therefore start the sequence with R- 2. 


The value of 


$y (%3) = 


(73) P (2) | 


may be obtained by direct calculation or trial 
and error and will yield both the minimum 
dy (#3) and the value of x, (a) at which it occurs. 
We may then proceed to find 4, (74) and wx, (x4) 
from (20) and so on. These could be tabulated 
as illustrated in an earlier paper [6]. However, 
it is interesting to observe that all this informa- 
tion can be included in a single graph such as 
Fig. 4. The curve marked I, is the curve of 
1/P (a) [continued at its extremities when the 
optimum temperature is outside (7',, 7*) by the 
curves 1/P* (x) and 1/P, (x)], the curve marked 
R has abcissa and ordinate 1/P (xp 
To show how this suflices, consider four reactors 
with feed x, and product 2, From as and the 
curve I, we know 1/P (23) (point A) and hence 
by drawing a horizontal line we have x, (point B). 
Again from x, and the curve I’, we have 1/P (2) 
(point C) and so x, (point D). Working up in this 
way from 2, to the curve I’, we have the complete 
optimum policy. The holding times of each stage 
can be quickly calculated from the area of the 
rectangles and the temperatures read off from the 
scale on the curve 1/P (2). 

Furthermore the curves I’, may themselves 
be found graphically by a simple construction 
rather than by a trial and error calculation of 


(20). Since the minimizing values of rp... . 
are in the open interval (x,.,.2,) they satisfy 

(21) 
Thus 


don 


by virtue of (21). To find the value of 2, which 
minimizes the right-hand side of (20) we have 


d | 
P (xp) + = 


P (xx) \ P 

P’ (rp) 
=P 


(22) 


Ms Xe Rs 2 


. 6. Construction for the optimum holding times 
of a sequence of stirred tanks. 


Consider the family of curves 


(2) 


P’ (2) 
h(a; %) P (2) (a) 


+¢ P2(2) (23) 
which may be drawn once and for all for given 
kinetics. Let this family be drawn for various ( 
in the same plane as the curve [,, 1/P (a). 


Chem. Engng. Sci. Vol. 13, Ne. 4. June, 1961. 


_ 


R. Aris 


Fic. 7. Construction for the curves J). 


(See Fig. 7, where four curves of the family are 
drawn). In case the temperature has to be held 
at r, or r* we interpret P (a) to mean P, or P*. 
When 2, is fixed the single-stage holding time ¢, 
is immediately determined for any 2,. The 


corresponding value x, (z7,) is always the given 
x, so that the curve I, is the horizontal line 
through A, the intersection of I, and the 
ordinate = Now for given a, the minimiz- 
ing 2, satisfies equation (22), ie. 


which means that the correct x, is the abcissa 
of the intersection of the curve h(a; ¢) for 
=a, with the point B. The point 25 = ¢ 
on the curve I’, has ordinate 1/K and so 
is obtained by dropping vertically from B to the 
point C on I’, and going horizontally from C to D 
on the ordinate 2 = (. Repeating this construc- 
tion with the various members of the family 
h(a, 0) for 0 (<a, a suflicient number of 
points to draw the curve I, can rapidly be found. 
I, is then constructed from I, in precisely the 
same way as I, was obtained from I. 

The value of this construction lies in the fact 
that I’, and the family h(x, 0) depend only on 


Fic. 8. The functions 1/P (y), hy (y) and A, (y). 
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the kinetics and the restriction on temperature ; 
from them such a diagram as Fig. 6 can be 
obtained for any c,. For the first-order reversible 
reaction we may use the variable y, as before, 
then 


1 
P(y) =y" (1 —y)? and h(y, = hy(y) — Che (y) 
where h, = (1 + p)y’ (1 — 


and hy =(p + y)y?*(1 — 


1/P (y), hy and h, are shown in Fig. 8 for the 
same values of p. Thus Figs. 2, 3 and 8 provide 
the material for a complete parametric study of 
the first-order reversible exothermic reaction. 

It may be remarked in conclusion that a 
diagram such as Fig. 6 could be constructed if 
the objective were to minimize some more com- 
plicated function F (@,, 0,...@,) than the sum 
However it is unlikely that simple 
graphical construction would exist and the 
functional equation corresponding to (20) would 
have to be solved by direct calculation. 

Since this paper was completed the work of 
F. Horn has come to my attention. The closed 
form of the optimal temperature and the resulting 
integrals was given by him in [7], where also he 
has tabulated some particularly useful cases of 
the integral. He has also used the same technique 
in [8] and [9]; for a more complete review of his 
work see [10]. 
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constant defined by (13) 

chemical species involved in reaction 

=1,...) 

extent of reaction 

extent of reaction in feed 

extent of reaction in reactor r 
concentration of A; 

concentration of A, in reference state 

tie 
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activation energies i = 1, 2 

minimum time to achieve conversion c, with 
feed cy in tubular or batch reactor 
minimum total holding time to achieve 
conversion c, with feed cp,, to R stirred 
tank reactors 

heat of reaction 
k,;* exp (— E,; RT), reaction rate constants 
R A cg, dimensionless maximum reaction 
rate 

E,/(E, — E,) 

maximum reaction rate 

r(c, T*), r(c, 

rate of reaction (dc dt) 

temperature 

upper and lower bounds of permitted 
temperatures 

time 

€1/ (eyo 

At 

stoichiometric coefficient in reaction 
a, A; 

exponents in rate expression 

pa; — B; 

Ad, 

holding time of stirred tank reactor r 
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R 
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constant defined by (15) 
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RT /(— AH) 


Chem, Engng. Sci. Vol. 18, No. 4. June, 1961. 


A 
A, 
(ey) 
tr (Cr+ 1) 
‘OL. AH 
13 k, 
96 P 
Pp 
Ric) 
R*, R, 
r(c, T) 
T 
T, 
t 
y 
a, 
By 
(2) 
(3) 
(4) 
205 


R. Aris 


(7] Horn F. Optimalprobleme bei kontinuerlichen chemischen Prozessen, Dissertation Techn. Hoscshcule Wien 
1958. 


{8} Horn F. Optimale Temperatur-und K trationsverliufe Second European Symposium on Chemical 
Reaction Engineering. Chem. Engng. Sci. 1961 14 77. 
[9] Horn F. and Trovrenter U. Chem. Ing. Tech. 32 1960 382. 


{10} Arts R. The optimal design of chemical reactors. New York. Academic Press. 1961. 


VOL. 
13 
96 


Chemical Engineering Science, 1961, Vol. 13, pp. 207 to 225, Pergamon Press Ltd., London. Printed in Great Britain 


New methods for the liquid phase continuous operation of 
thermogravitational columns with reservoirs 
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Abstract—-Experimental and theoretical results are presented for four new methods for the 
liquid phase continuous operation of thermogravitational columns with reservoirs. Two of the 
methods for continuous operation provide techniques for experimental study and refinement of 
thermogravitational column theory. The developed theory for the four new methods is in quali- 
tative agreement with the experimental results obtained for the binary liquid system, n-heptane 
benzene. 


Résumé — Présentation des résultats théoriques et expérimentaux de quatre nouvelles méthodes 
pour la marche continue en phase liquide de colonnes A circulation thermique avec réservoirs, 
Deux de ces méthodes apportent des données techniques pour l'étude expérimentale et le per- 
fectionnement de la théorie des colonnes 4 circulation thermique. La théorie des quatre 
nouvelles méthodes est en accord qualitatif avec les résultats expérimentaux obtenus pour le 
systéeme liquide binaire n-heptane-benzéne. 


Zusammenfassung — Experimentelle und theoretische Ergebnisse werden fiir 4 neue Methoden 
des Betricbs von kontinuierlichen Flissigphase Thermogravitationssiulen mit Reservoirs mit- 
geteilt. Zwei dieser Verfahren fiir kontinuierlichen Betrieb liefern Hinweise iiber experimentier- 
technische Studien und fiir eine Verfeinerung der Theorie der Thermogravitationssiulen. Die 
entwickelte Theorie fiir die vier Verfahren stimmt qualitativ mit den experimentellen Ergebnissen 


tiberein, die fiir das fliissige binidire System n-Heptan-Benzol erhalten wurden. 


CONSIDERABLE interest has been shown in the 
thermogravitational column with the result that 
quite an extensive body of literature has developed 
on the subject, as is evidenced by several reviews 
{1 -38]. The majority of the experimental work 
on the thermogravitational column has involved 
operation of the column in what is termed a 
“batch” manner. Operation of the thermo- 
gravitational column in a “* continuous manner 
has been the subject of more recent investigations. 

Batch” operation of a thermogravitational 
column consists of charging the column initially 
with the mixture to be separated and permitting 
no additional material to cross the boundaries of 
the column during the course of the separation. 
The theory for batch operation has been treated 
by many investigators [4-9], resulting in at 
least a correct qualitative description of experi- 
mental results. 


There have been a number of relatively recent 
experimental studies of the continuous operation 
of thermogravitational columns [2, 3, 10— 12}. 
The increasing industrial importance of separation 
processes has resulted in a substantial patent 
literature relating to design modifications for 
continuous process separation equipment utilizing 
thermal diffusion e.g. [13-27]. Most of these 
patents represent some alteration of thermo- 
gravitational column construction or operation 
with the view of increasing separation efliciency. 
The theory for the operation of a “* centre-fed ” 
thermogravitational column has been considered 
by several investigators [2, 3, 7, 10]; the term 
centre-fed describes the continuous operation 
in which feed is introduced at some vertical point 
of the column (usually the vertical centre, thus 
“centre-fed’) and product streams are with- 
drawn from the top and bottom of the column 


. 
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DescRIPTION OF THE NEW METHODS FOR THE 
LIQUID PHASE CONTINUOUS OPERATION OF A 
THERMOGRAVITATIONAL COLUMN WITH 
RESERVOIRS 


Experimental studies, excluding patent work, 
have considered only one continuous method for 
operating a thermogravitational column, “ centre- 
fed” operation. It is the purpose of this paper 
to describe four new methods proposed by the 
author for continuous operation of thermo- 
gravitational columns. While the new methods of 
operation in altered form might be applied to 
gases, only the operation with liquid mixtures is 
considered here. 

The apparatus used in the present study is 
described as a parallel plane plate thermogravita- 
tional column with reservoirs at the extremities 
of the vertical dimension. The reservoirs have 
volumes which are effectively large with respect 


Fic. 1. A flat plate thermogravitational column with 
reservoirs : co-ordinate system and column dimensions 
Legend : 


4 = Column heat transfer area height 
L_ = Column length 
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to the volume contained between the parallel 
plates ; the liquid filling the reservoirs communi- 
cates freely with the liquid between the parallel 
plates (see Fig. 1). 

Thermogravitational columns with reservoirs 
have been extensively employed in experimental 
studies of the unsteady-state batch operation of 
thermogravitational columns. The presence of 
reservoirs simplifies the mathematical description 
of the unsteady-state problem and has a practical 
effect in making possible greater accuracy in 
experimental work. 

This study is concerned with the liquid phase 
continuous operation of a flat plate thermo- 
gravitational column with reservoirs. As an aid 
in understanding the description of the experi- 
mental and theoretical work which follows, each 
of the four new methods of continuous operation 


will be described brietly. 


Two-feed cocurrent operation 

The two-feed cocurrent operation of a flat 
plate thermogravitational column with reservoirs 
consists of the introduction at the same horizontal 
end of the column of liquid feed, not necessarily 
of the same concentration, at the upper and lower 
reservoirs. Liquid is withdrawn from the upper 
and lower reservoirs at the other horizontal end 
of the column. The flow pattern is indicated in 
Fig. 2 (a). 

Introduction and withdrawal rates of liquid at 
the upper reservoir are adjusted so that there is 
a constant volumetric flow rate in the upper 
reservoir. Similarly a constant volumetric flow 
is maintained in the lower reservoir by proper 
adjustment of introduction and withdrawal rates. 

In a thermogravitational column designed for 
operation with liquids the perpendicular distance 
between the hot and cold walls is usually in the 
range 10°? 10-' cm. In this study the spacing 
between hot and cold walls was 0-0376 cm. The 
reservoirs were each semi-circular in section with 
a diameter of 0-952 cm. The significance of these 
dimensional values is that the pressure drop for 
flow between the parallel plates at isothermal 
conditions is of the order of a hundred times that 
for flow through the reservoirs. 

Owing to this pressure-drop relation, the only 
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flow occurring during two feed cocurrent operation 
between the parallel plates is the natural convec- 
tion arising from the temperature gradient ; 
forced convection occurs only in the reservoirs. 
There are then two separate transverse flows in 
the column: natural convection in the + y 
direction arising from the temperature gradient 
and forced convection in the +2 direction 
(cf. Fig. 2a) in each of the reservoirs arising from 
an externally imposed pressure gradient. 

Calculations and experimental flow studies with 
the apparatus [28] showed that the flow conditions 
described above are substantially correct. Agree- 
ment between experimental and theoretical results 
lends further support for the described flow system. 

There are several patents on cocurrent-type 
operation [15, 17, 29). The various patented 
methods differ from that of the present study in 
that pesitive flow of material in the present 
method occurs only in the reservoirs. Some of 
the previous methods [15, 17] utilized plane plate 
construction without reservoirs ; liquid was forced 
between the parallel plates with a_ resultant 
interference between foreed convection in the 
+ z direction and natural convection in the + y 
direction. 


Two-feed countercurrent operation 


This method of operation uses the same flow 
principles described above for two-feed cocurrent 
operation. The difference in two-feed counter- 
current operation is that the feed streams to the 
two reservoirs enter at opposite ends of the ap- 
paratus as shown in Fig. 2(b). The product 
streams are also withdrawn from opposite ends of 
the column so that the flows in the two reservoirs 
move in a countercurrent manner to one another. 

Just as for two-feed cocurrent operation 
natural convection occurs between the parallel 
plates and forced convection occurs in the 
reservoirs. Separate, transverse flows occur in 
the apparatus for two-feed countercurrent opera- 
tion. 


Single-feed cocurrent operation 

Single-feed cocurrent operation denotes a 
method of operation in which liquid mixture is 
fed into only one of the two reservoirs. This 
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Fic. 2. Schematic flow diagram of the four methods 


for the liquid phase continouus operation of a thermo- 
gravitational column with reservoirs. 
Legend : 

w = Mass fraction of that component of a binary 
mixture migrating to “* hot” wall and concen- 
trating toward the top of the column. 

s = Product mass flow rate. 


procedure is illustrated in Fig. 2(¢) where liquid 
is shown fed to the upper reservoir and product 
withdrawn from each reservoir at the other 
horizontal end of the column. 

The pressure drop relation noted above in 
the discussion of two-feed cocurrent operation 
is of importance here. When liquid is withdrawn 
from the upper and lower reservoirs as shown in 
Fig. 2(c) the total product flow from the lower 
reservoir has come from the upper reservoir, and 
this flow was distributed uniformly over the 
horizontal length (z direction) of the column. 
Thus the volumetric flow rate determined from 
point to point between the parallel plates at 
isothermal conditions would have a constant value 
for fixed feed and withdrawal rates. 

The difference between two-feed operation and 
single-feed operation now appears. With single- 
feed operation there is no longer a separation of 
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natural and forced convection, Between the par- 
allel plates there occur both natural convection 
(in the + y direction) and forced convection (in 
either, but not both, the + or — y direction). In 
addition with single feed operation the volumetric 
flow rate in the reservoirs is not a constant but 
is, to an approximation, a linear function of the 
horizontal dimension, =. 


Single-feed countercurrent operation 


Single feed countercurrent operation is similar 
to the single-feed cocurrent method described 
above, except that here liquid is withdrawn at 
opposite ends of the column. Fig. 2(d) illustrates 
the situation in which feed is introduced into the 
upper reservoir at one horizontal end (z — 0) 
and liquid is withdrawn from the upper reservoir 
at the other horizontal end of the column (2 — L) 
and from the lower reservoir at the horizontal end 
where the feed was introduced (2 0). 


EXPERIMENTAL APPARATUS AND PROCEDURE 


Only brief discussion will be given the 
experimental apparatus and procedure; full 
details can be found elsewhere [28]. 


Experimental apparatus 

The flat plate thermogravitational column with reservoirs. 
The column was constructed of two rectangular pieces of 
Alcoa tool and jig plate, both plates having dimensions of 
approximately 34in. I2in. gin. The plates had been 
machine polished. The edges of the two plates were 
milled and one groove, semicircular in section and having 
a diameter of jin. was machined on the surface of each 
plate (see Figs. 3a and 4). 

When placed in operating alignment the groove in the 
“hot” plate was at the bottom of the column and the 
groove in the “ cold” plate was at the top of the column. 
The parallel surfaces between the two grooves had a height 
of 19-33 em. 

Sealing and formation of the space between the two 
plates was accomplished with a gasket cut from a sheet 
of Teflon, 34in. « 12in. « 0-0150in. | 0-0005in. as shown 
in Fig. 3(b). Thin vertical strips at the positions indicated 
in the Fig. were found necessary for the achievement of 
a suitable plate spacing. For the column dimensions and 
the methods of operation studied these strips had a 
negligible influence on the experimental results. 

In assembling the column the gasket was sandwiched 
between the surfaces of the two flat plates, the cooling 
conduit was placed on the back of the * cold ~ plate and 
the whole assembly was held together with sixty-five 
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Fic. 4. Details of construction of thermogravitational 
column with reservoirs, 


Detail section of upper reservoir 


Legend : 
S.T. — Sample tap 
FP = Feed or product access 
Gs = Gasket 
Seale: lin, — 2in. 


24in. C-clamps placed at determined positions, Tests under 
pressure revealed no leaks in the assembly. 

A satisfactory plate spacing was achieved by a system 
of pressure bolts which could be tightened against the back 
of the cold plate, thus pushing the cold plate toward the 
hot plate and eliminating a bowing tendency found when 
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the plates were clamped together. Three small 
(yin. in. a) aluminium spacers were used together 
with the Teflon strips of the gasket to maintain a given 
plate spacing ; “a” is the plate spacing. 

Plate spacing was determined by injecting a modified 
liquid plaster of Paris mixture into the assembled column 
and allowing the mixture to harden. Through a trial-and- 
error procedure the correct torque wrench setting for the 
pressure bolts was determined to give a desired plate 
spacing. In a typical plate spacing determination, out of 
26 measurements spaced over the heat transfer surface, 
20 were within + 0-0002in. of 0-0152in. ; the remainder 
deviated at most by + 0-0005in. 

In order that all experimental runs would be self-con- 
sistent the column was assembled in the manner described 
above and sufficient time was allowed for cold flow of 
materials to take place (judged as 2 months). All ex- 
perimental runs were then conducted with the assembly 
with checks on plate spacing furnished periodically by batch 
steady-state separation determinations. After all runs 
were completed a plaster of Paris plate spacing determina- 
tion was made. For the experimental runs reported here 
the periodic steady state separation determinations gave 
the same value over the time of the runs and the final 
plate spacing determination based on 26 measurements 
gave a viaue of 0-0148in. | 0-0005in. This is the plate 
spacing value used in all subsequent calculations. 

Cooling system, The “cold” wall temperature was 
maintained by flowing water through a conduit which 
fitted over the back of the “cold” plate. The area 
perpendicular to the flow direction was 9-5in. « fin. 

Provision was made for mixing water available at 
80 °C with water available at a minimum of 13 °C so that 
the * cold * wall temperature could be adjusted to a desired 
value, 

The cooling water was passed through a series of two 
surge tanks to eliminate small pressure fluctuations 
produced by the control valves. Connection was made 
to the cooling conduit by two jin. inlets ; cooling water 
exited from the conduit through two jin. outlets. Vibra- 
tion was minimized by the use of rubber hose connections. 

During operation the “* cold * wall temperature over the 
plate surface deviated from an average value by + 0-5 °C, 

Temperature of the “cold” wall was determined by 
six iron-constantan thermocouples distributed over the 
“cold” plate. The junctions were fitted into holes 
extending to within gin. or less of the “cold” wall. 
The thermocouple leads were connected to a Leeds and 
Northrup Speedomax recorder. Recorder and thermo- 
couples were calibrated against a Leeds and Northrup 
Precision potentiometer. 

Heating system. The “hot ™ wall temperature was 
maintained by a specially constructed electrical radiant 
heater. By varying the radiant properties of the back of 
the “ hot " plate, end losses and convection effects could 
be corrected. 

A General Radio Corp. Type 50B, 230 V, 7 kVA, 50-60 


cycle variable reactance provided a variable voltage to 
the heating unit. 

Using the heating system just described it was found 
that the “ hot” wall temperature could be maintained 
to within + 0-7 °C of a given value over the surface in the 
range 49-59 °C, The “hot” wall temperatures were 
measured by six iron-constantan thermocouples placed 
with junctions jin. or less from the “ hot” wall. The 
thermocouple leads were connected to the temperature 
recorder described above. 

Feed and product system. The flow system used in 
this investigation is shown schematically in Fig. 5. Feed 
mixture was pumped by a Milton-Roi “* miniPump ” 
from a stirred collecting vessel to a constant-head feed 
vessel; overflow from the feed vessel returned to the 
collecting vessel. 

Liquid feed was conducted from the feed vessel through 
two feed rotameters (for two-feed operation) with regula- 
tion supplied by two feed needle valves. Product from 
the thermogravitational column was passed through 
two product rotameters with regulation supplied by two 
product needle valves and 3in. lengths of copper tubing 
packed with 200 mesh glass spheres. 

Feed and product rotameters were obtained from the 


Flow system for flat plate thermogravita- 
tional column with reservoirs. 


Legend : 


Feed valves 

Feed rotameters 

Thermogravitational column with reservoirs 
Product valves 

Product rotameters 

Packed resistance tubes 

Product mass flow rate measurement 
Collecting tank 

Magnetic stirrer 

Feed pump 

Constant head feed tank 

Overflow return line 
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Manostat Corp., Cat. # 9143b. Modification of the tubes 
to eliminate parallax made it possible to reproduce a 
given flow rate in the rotameters to + 0-1 per cent or 
better. 

The liquid mixture studied in this investigation contacted 
only the following materials: glass, stainless steel, brass, 
copper, tin, aluminium, Teflon and sapphire. 


Experimental procedure 

Feed mixture. The binary mixture, n-heptane-benzene 
was chosen for this study. 

Benzene (Reagent grade) was obtained from the Fisher 
Scientific Co. and had a measured refractive index at 
25 °C of 1-4979 which agrees with the quoted literature 
value [30]. n-Heptane was obtained from the same vendor 
and had a refractive index at 25 °C of 1-3852 which agrees 
with the quoted literature value [30]. No further purifica- 
tion of the chemicals was made. 

Flow rate measurement. A timed sample of each product 
flow was taken in a glass weighing bottle fitted so as to 
prevent evaporation. The timing of the collected sample 
was performed with an electric timer from the Standard 
Electric Time Co. The collected, timed samples were 
weighed on a Christian Becker Chainomatic analytical 
balance. 

Composition measurement. Composition of the liquid feed 
and product streams was determined by measurement of 
the refractive index at 25 °C with a Bausch and Lomb 
Abbe Refractometer. 


THEORETICAL DESCRIPTION 


Description of the four methods of operation 
in terms of experimental quantities depends on 
solution of the equations of conservation of mass, 
momentum and energy for the physical system. 
This solution is accomplished here by development 
of the * transport * equation through a procedure 
originated by Furry et al. [7] for gases and ex- 
tended to liquids by pr Groor [8]. The * trans- 
port’ equation is a statement in differential 
form of the net rate at which a particular com- 
ponent of a binary mixture is being “ transported” 
in the vertical direction in a thermogravitational 
column. 

Other methods have been developed for the 
solution of the conservation equations in describ- 
ing the operation of the thermogravitational 
column. The various methods of solution have 
been discussed in a comprehensive review by 
Lonemire [2]. 
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Derivation of the transport equation for the liquid 
phase continuous operation of a flat-plate thermo- 
gravitational column with reservoirs 


The equations of conservation of mass, momen- 
tum and energy have been discussed by Birp [31] 
whose nomenclature will be used in this dis- 
cussion. 

Fig. 1 presents the description of the geometry 
of the flat plate thermogravitational column 
with reservoirs. The right-handed Cartesian 
co-ordinate system has its origin at the lower 
edge of the “ hot ” wall. The + z axis runs along 
the lower edge of the * hot * wall. “ L” denotes 
the length of the heat transfer area along the 
2 axis. “a” is the distance along the +2 axis 
to the “cold” wall. “A” is the height of the 
heat transfer area along the + y axis. 

The following assumptions are made in the 
development which follows : 


1. Steady state 
2. Binary mixture 
3. No chemical reaction 


Laminar, Newtonian, flow 
v= y(t); haz p= ply); 
O< y< h,a,z. 


v=v,; h<y2,% 2%. 
p const; 
O>y h< yw, 
5. Only thermal and ordinary diffusion contri- 
bute to the mass flux. 


6. Thermal conduction is the only important 
mode of heat transport 
T=T(2r); O<y< h,z,2. 
T =(T O>yh< y, 2, % 


where v denotes the velocity, v, is the average 
velocity in the upper reservoir, »; the average 
velocity in the lower reservoir, p is the pressure, 
T the temperature, and (7',),, the average tem- 
perature in the reservoirs. 

The transport equation will now be developed 
for two-feed operation ; assumptions described 
later permit a simple modification of this develop- 
ment so that single-feed operation is described. 
With the assumptions noted above the solution 
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of the energy balance equation for the region 
0< y.< is 

dT Ty AT 

dr a 

The solution of the equation of motion becomes 

for the same region with the usual assumptions 
and boundary conditions : 

(cos 0) Bp g AT 


vy (2) = - — a); 


(1) 


O<y<h,a,z (2) 


a third-degree equation for the velocity. 

For the assumptions above the mass flux 
expression for component “1” of the binary 
mixture is 


where j, is the mass flux vector referred to a co- 
ordinate system moving with the mass average 
velocity, [g (em* sec~!)], p is the mass density of 
the solution, (g/em*), D,, is the mutual diffusion 
coeflicient for a binary mixture, (cm? sec), w, the 
mass fraction of that component of the binary 
mixture travelling to the “hot” wall and en- 
riched at the top of the column (component “ 1”) 
a, is the Soret coefficient of “1”, (degree~*). 

Inserting equation (6) in the equation of 
continuity gives with appropriate assumptions 
for the region 0 < y < h, x, 2. 


Mw, 
dy? 
AT) 
where it has been assumed that the dependence 
of p, Dy, and o, on x, y and = can be accounted 


by assuming average values. A solution of equa- 
tion (7) is sought for the boundary conditions 


+ = 0 a=0,a (5) 


Conditions for the y and z dimensions will be 
introduced below. 

In developing the “transport equations 
Furry et al. [7] and pe Groor [8] neglect for the 


moment the diffusion in the y direction and later 
in the development account for this term. It 
will be found that only for very small values of 
the plate spacing, a, does this neglected term 
contribute to the * transport.” 

Strictly the term accounting for diffusion in 
the z direction cannot be neglected at this point 
unless it is negligible. Here diffusion in the z 
direction will be assumed negligible on the same 
basis that diffusion in the y direction can be 
neglected ; this basis is that the transfer of “1” 
by the convective flows is much greater than the 
transfer of “1” by ordinary diffusion. 

These diffusion terms may be introduced by 
writing for the net rate of transport of component 
“1” up the column, m (2, y, z), at any point, z 

a 


dm =) [re y, 2) dx 
0 


where + (7, y, 2) is the concentration diffusion flux 
in the z direction at a point in the column in the 
region O< Since the change in 
concentration in the 2 direction is small, an 
average value, 7(y, 2), for this diffusion flux can 


be assumed 


dm (x, y, 2) 
= ?(y,2 -pD 
dy (y, 2) P 12 
where @, denotes the mass fraction of “1” 
averaged over x. It follows that 


m (2, y, %) = 0, 2) + 


0 


This equation indicates that if concentration 
diffusion in the z direction is negligible the net 
rate of transport of component * 1” is a function 
only of 2 and z. By a similar procedure diffusion 
in the y direction can be accounted by including 
in the net rate of transport of “1” up the column 
the transfer by diffusion in the y direction. 
Assume that concentration diffusion in the z 
direction is negligible compared with the rate at 
which “1” is transported by free convection. 
With this assumption one then obtains the result 
of Furry et al. [7] and pe Groor [8] which states 
that in the 7-y plane at any point z the net rate 
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at which component “|” is moving up the 


column, m (z) is: 
m (z) = 


_ [1 (Bpg\? pa? 
(524) costo + x 


for the region 0 < y < h,a, 2. The usual assump- 
tions have been made in obtaining this equation. 

The difference between equation (6) and the 
“transport” equations of Furry ef al. and 
pr Groor is that equation (6) defines the * trans- 
port” in a plane perpendicular to the = axis in 
a system in which = is an independent variable ; 
[7] and [8] did not consider z as a variable in their 
developments. 

Introducing the notatior of the “ transport 
constants ” H* and A* into equation (6) 


dw, 
oy 
where H* and A* are defined by equation (6), 
Equation (7) is the form of the “ transport ” 


equation to be used in the developments for two- 
feed cocurrent and countercurrent operation, 


(7) 


m (2) -H* 


Separation equations for two-feed cocurrent and 
countercurrent operation 


Two-feed cocurrent operation, This method of 
operation has been described previously and is 
illustrated schematically in Fig. 2(a). 

In what follows it is assumed that the same 
degree of mixing prevails along the entire length 
of the reservoirs. This assumption applies to the 
other methods of operation to be discussed. The 
assumption holds since temperature conditions 
are the same over the length of the reservoir and 
there is a relatively large natural convective 
velocity in the reservoir which describes, in column 
operation, a helical path about the 2 axis. 

In the development of equation (7) it has been 
assumed that there are separate, transverse flows 
in the column corresponding to the forced convec- 
tion in the reservoirs and natural convection 
between the parallel plates. With this assumption 
a material balance can be written about any 


element Az of the reservoir length, L. For two- 
feed cocurrent operation, the reservoir material 


balances are 
(8) 


(9) 


where average values for the mass densities of the 
solutions in the reservoirs have been assumed to 
account for the variation of p with 2. s, and s, 
are the average mass flow rates in the upper and 
lower reservoirs respectively (g/sec) and ,, and 
w,, are the mass fractions of component “17° in 
the upper and lower reservoirs. 

Since the “ transport ” m (z), at any = is inde- 
pendent of y (if concentration diffusion in the z 
direction is negligible) the integrations of equation 
(7) over y and z= may be performed separately. 

Equation (7) can be integrated over the y 
direction as it is written. However, the non- 
linearity does not permit the following integration 
in the z direction in closed form. The problem of 
the integration of the “ transport " equation has 
received attention from Jones and Furry [82] 
who recognize two approximations, which for the 
present purpose will give expressions for the separ- 
ation over the column in a convenient form. 

The following approximations are made : 


Wy Wy, O38 < Ww, <= 07 
4; OF (10) 


for concentrated solutions where ¥, and #, are 
average mole fractions of “1° and “2. For 
dilute solutions the linear approximation 


- aa +b (11) 


will be used. For all dilute solution experimental 
runs of this investigation the constants are taken 
to be a’ = 0-7, b = 002. 
Equation (7) is integrated over y with the 
boundary conditions 
yooh (12) 


Equations (8) and (9) and the integrated form 
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of equation (7) are solved for the boundary 
condition : 
A=2,, —2#,4=90, z:=0 (13) 
The solution is for the approximation of equation 
(10) 
— = 
= H*h 


For the dilute linear approximation equation (11), 
the solution is 


Oy = + 5) x 


li 1 — exp (— 4,) 
1 +(8,/8;) exp(— 3) 


_ [1 + (s,/s;) exp (— ¢,)]= | 
(s, a’H*) (1 — exp(— 


x (1 + “4 (15) 


where: ¢, a@H*h/K* (16) and w, denotes the 
mass fraction of * 1” in the feed mixture. 
Two-feed countercurrent operation. The develop- 
ment of separation equations for two-feed 
countercurrent operation parallels the develop- 
ment for two-feed cocurrent operation, 
The reservoir material balances are 


« exp | 


(17) 


(18) 


Equations (17), (18) and (7) are to be solved for 
the boundary conditions 


(19) 


The separation equation for the dilute approx- 
imation, equation (11), is 


[wp + (b/a’)) [1 — exp (¢)) 
(8,/s)) exp (— BL) — exp (¢,) 
x [1 —exp(— Bz)] (20) 


— Op = 


where ¢, is defined in equation (16), Also 


[= (4,) — (8,/8)) 
8, exp (¢,) — 1 


Separation equations for single-feed operation 


(21) 


In single-feed operation there is a positive flow 
of fluid in the y direction produced by a pressure 
difference between the two reservoirs. Thus, in 
the space between the parallel plates forced con- 
vection is superimposed on the free convection 
arising from the temperature gradient. 

For a unidirectional velocity field produced by 
a very small additional pressure difference 
Jones and Furry [82] point out that this field may 
be superimposed on the free convection field. 

If the further assumption is made of a small 
variation of the concentration of “1” in the y 
direction (this assumption was implied in the 
derivation of the “transport” equation) the 
only change in the “ transport ” of “1” up the 
column for single feed operation takes the form 
of an additional term to account for the positive 
transport of “1” in the + or — y direction 
caused by forced convection. 

It was pointed out that to an approximation 
the total flow out of the feed reservoir to the 
other reservoir is distributed uniformly over the 
total column length, L. With the assumptions 
discussed above the “ transport ’’ equation for 
single-feed operation becomes 


m'(z) = H* + K* +r@, (22) 


where m’ is now the gem~! sec~! of species * 1” 
transported down the column; r is the positive 
total flow down the column (g sec?). 

Separation equations for single-feed cocurrent 
operation. The development of separation equa- 
tions for single-feed operation is analogous to the 
procedure followed for two-feed operation. 

The reservoir material balances have the form 


—m' (23) 


(24) 


where s, is the total feed rate to the column 
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(g/sec). The solution to equations (22), (23) and 
(24) is subject to the boundary conditions on y. 


(25) 


my, Yoh 


== @ (26) 


In order to obtain equations in closed form the 


following assumptions are made in solving 
equations (22), (23) and (24) : 

27 


for the case of concentrated solutions and 


op (a’ r) ‘| at 


28 
K (28) 


for dilute solutions. 
The separation equation for the concentrated 


solution approximation, equations (10) and (27), 
is 


So = In (1 (29) 


(4) 


For the dilute linear approximation, equations 
(11) and (28), 


VWs 
“me 


the separation equation for dilute solutions is 
obtained, 

Separation equations for single-feed countercurrent 
operation. The development of separation equa- 
tions for this method of operation follows the 
procedure given for single-feed cocurrent opera- 
tion. 


The reservoir material balances are 
da, 
(% — rz) ~ Torn, m 
r(L 2) Sow Tay = — 


When the assumptions of equations (27) and (28) 
are made the developed separation equations for 
this method of operation are identical to equations 
(29) and (30). 
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EXPERIMENTAL RESULTS 

Experimental runs were made for the four 
methods of operation using the binary system, 
n-heptane—benzene. Additional experimental data 
were taken for the unsteady state batch operation 
of the column of this investigation ; these data 
are reported elsewhere [28] and are not given 
here. 


Two-feed cocurrent operation results 

The most direct test of the various assumptions 
given above in the discussion of the theory of 
two-feed cocurrent operation is provided by a 
comparison of the experimental data for two-feed 
cocurrent operation with the corresponding sep- 
aration equations, equations (14) and (15). 

For the concentrated solutions the separation 
equation, equation (14), may be arranged to give 
the following indicated relationship : 


A 
Cin(t x) (31) 
where C, is a constant : 
C, (h/ (32) 


and s, and s,; are the mass flow rates in the upper 
and lower reservoirs, g sec~', A is the difference 
in concentration of species 1” in the upper and 
lower reservoirs for given mass flow rates and A, 
is that steady state separation determined by 
extrapolation of experimental data for two-feed 
cocurrent operation to zero mass flow rate, and 
as also determined from experimental batch 
separation data. 

For the case of dilute solutions the separation 
equation, equation (15), may be arranged ; 


where C, is a constant : 


[1 — exp (— 4,)] 
a H* + exp(— 


and $ is the average of the flow rates for “ sym- 
metric * runs in the upper and lower reservoirs. 
The term “symmetric” indicates that the 
experimental run was performed with approx- 
imately equal mass flow rates in the upper and 


(33) 


(34) 


VOL. 


New methods for the liquid phase continuous operation of thermogravitational columns with reservoirs 


Eqni4, log ( 4) 


Fic, 6. Agreement of experimental data with linear 
relation predicted by equation (14) for two-feed co- 
current operation. 


Legend : 


tp 0-500 0-500 


Runs TF-SCO-K 
1, 3,4 15 
a5 6, 7 
A =4, — separation 
Ay = Separation, zero flow rate 
a2 = Mole fraction n-heptane in reservoir product 
s — Product mass flow rate 


lower reservoirs, — a situation realized to + 0-1 per 
cent of the average for the larger flow rates. 
Figures 6 and 7 present the experimental sub- 
stantiation of the indicated relationships of 
equations (31) and (33). Fig. 6 presents experi- 
mental data for a feed concentration of 0.500 
mole fraction n-heptane (termed the K series) 
and the indicated three different temperature 
differences. Fig. 7 shows experimental data for a 
feed concentration of 0-150 weight fraction 
n-heptane (termed the * D” series) at a temper- 
ature difference of 31-6 °C. It can be seen that 
the agreement with the indicated relationship 
is good, At lower reservoir mass flow rates the 
experimental data do not follow these straight 
lines, since for these low flow rates the separation 
becomes large and therefore the linear approxima- 
tions, equations (10) and (11), and the assumption 


Eqn 


Fic. 7. Agreement of experimental data with linear 
relation predicted by equation (15) for two-feed 
cocurrent operation. 
Legend : 
Two-feed cocurrent experimental point : 
Runs TF-SCO-D, 1-4 
: 316°C 
: 
= Steady state separation, zero flow rate 
= Ws — w;, Separation 
= Mass fraction n-heptane in reservoir product 
= Arithmetic average of upper and lower product 
flow rates = $(s, + s,) 


of constant mass flow rate in the reservoirs fail 
to describe the physical situation ; also restric- 
tions applied in the derivation of the * transport ” 
equation, equation (7), are violated. 

The constant in equation (31), the slopes of the 
straight lines of Fig. 6, contains the column height, 
h, length, L and the “ transport constant,” A*. 
Therefore K* can be measured by direct experi- 
mentation. This is analogous to the procedure 
which has been used [3, 33] for determination of 
transport constants from unsteady-state batch 
operation data, H* can be calculated in turn from 


(35) 


using the experimentally determined zero flow 
rate steady state separation Ay. Another method 
for determination of H*, which in general will not 
agree with the value determined from equation 
(35), since the linear approximation is introduced 
in deriving equation (35), is to calculate H* from 
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Table 1. 


Comparison of experimental and calculated values of * Transport Constants".” 


K* cale. 
equation (6) 
44 
+ x 


10-4 
1-4 


318 | 
| 


AT 
T = 307-3 °K 


TF-SCO-K 
AT € 
T = 304-9 °K 


TF-SCO_K 
AT = 395°C 
T = 3124 


47 x 


TF-SCO_D | 44 «10-4 
+02 x 


AT = 316°C x 104 
T = 3061 K | 


— 183 x 10-5 
+ 0038 


— 677 lo 
+ 006 « 10-5 


‘. K series transport constants calculated from mole fractions ; D series transport constants calculated from 


mass fractions, 


2. Code explanation: K, Feed concentration 0-500 mole fraction n-heptane in benzene ; D, Feed concentration 
0-150 mass fraction n-heptane in benzene ; TF, Two feed operation ; CO, Cocurrent operation ; S. Symmetrical operation, 


approximately equal product mass flow rates. 


*, H* values for K series calculated using equation (35); H* values for D series caleualted using equation (15). 


the solution of the “transport” equation, at 
zero flow rate and steady state. The solution is: 


ayy |- he 


Wis | (36) 


For consistency the calculation of 1/* by equation 
(35) is to be preferred. 

By a similar procedure for the dilute solutions 
H* and K* can be calculated from the slope of 
the line of Fig. 7 and equation (33). 

It is well known that experimentally determined 
“transport constants” do not agree in general 
with values calculated from the defining equation, 
equation (6) for example. The disagreement is 
particularly serious for liquids and has not yet 
been satisfactorily explained. 

Table 1 is a summary of the comparison of 
calculated (equation 6) and experimental values 


of the “ transport constants.” Physical property 
data for the system n-heptane-benzene used in 
calculating H* and A* from equation (6) are from 
Refs. [10] and [34]. In agreement with com- 
parisons for other liquid systems the Soret 
coetlicient determined in the thermogravitational 
column is much larger than the coefficient 
determined by the methods in which all convec- 
tion effects are eliminated, For the present 
system the comparison of the Soret coeflicient 
determined by the thermegravitational column 
with values determined by the unsteady-state 
cell method [34] is made at two difierent con- 
centrations in Table 2. 

In further consideration of the other operating 
methods of this investigation the experimentally 
determined “ transport constants” of Table 1 
will be used in the various equations which have 
been developed. 
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Table 2. Comparison of Soret coefficients deter- 

mined from thermogravitational column measure- 

ment in this investigation with Soret coefficients 

determined by measurement of pure Soret effect. 
System: n-heptane—benzene 


Average temperature : 307-2 °K 
Pressure : 1 atm 


Uss Gl 
measured by 
et al. [34) 


Mass fraction col. this 
n-heptane | investigation! 


0-562 49 — 1-37 4 0-04 (10-4) 
0-150 714 03(10-%) — 302 4 0-04 (10-%) 
| 
1 Sign is negative if n-heptane goes to hot wall. 
2 Calculated from mass fractions. 


Figures 8 and 9 represent all the data for two- 
feed cocurrent “symmetrical” operation. For both 
the concentrated solution (0-500 mole fraction 
n-heptane) and the dilute solution (0-150 mass 


seca 


Fic. 8. Upper and lower reservoir product concentra- 
tion difference vs. sum reciprocal product rates : 
two-feed cocurrent operation. 


Legend : 


4 | | oO 


AT(C) 39-5 31-8 28-5 
Runs TF-SCO_K | TF-SCO-K | TF-SCO_-K 


Eqn. 


a = Mole fraction n-heptane in product 
s = Reservoir product mass flow rate 


fraction n-heptane) it can be seen that equations 
(14) and (15) predict larger separations at the 
lower flow rates than found experimentally. The 
assumptions made in deriving the separation 
equations are valid only for small separations ; 
the disagreement at the lower flow rates therefore 
is to be expected. 


Fic. 9 Upper and lower reservoir product concentra- 
tion difference vs. reciprocal average product rate : 
two-feed cocurrent operation. 
Legend : 
© = Runs TF-SCO_D, 1-7 
equation (15): @ =0O7, b= 0-02 
AT = 316°C 
wp = 0-150 
w = Mass fraction n-heptane in product 
¥ = Arithmetic average of reservoir product mass 
flow rates = (8, + 8;) 


Figures 10and 11 present the experimental results 
of * asymmetric ” (unequal flow rates in the two 
reservoirs) two-feed cocurrent operation. The 
experimental results are indicated by the solid 
curves which represent separations obtained at a 
fixed lower reservoir product mass flow rate while 
varying the upper reservoir product mass flow 
rate. Comparison is made between the experi- 
mental curves and those predicted by equations 
(14) and (15). As expected, the agreement becomes 
poorer with increasing separation. 

It is of interest to compare the experimental 
separations found with two-feed cocurrent opera- 
tion with those found in “ centre-fed * operation 
of a thermogravitational column. Fig. 12 presents 
such a comparison. The two upper curves repre- 
sent the experimental curves obtained in the 
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Fic. 10. Feed and lower reservoir product concentra- 
tion difference vs. upper reservoir flow rate: two-feed 


cocurrent operation. 
Legend : 
AT = 317°C 
Runs TF-USCO_K, 1-18 
8, 0001645 g sec 
8; ~ 0-003134 g/see 
8; = 0-005649 g sec 
~— equation (14) 
« = Mole fraction n-heptane 
s = Product mass flow rate 


4. 


0-08 O10 


Feed and lower reservoir product concentra- 


tion difference vs. upper reservoir flow rate : two-feed 


cocurrent operation. 
Legends : 

AT 316°C 
Runs TF-USCO-D 1-10 
8; = 0001683 g sec 
8; = 0-002541 g/sec 
© 8; = 0-004664 g/sec 
equation (15) a = 0-7; b= 002 
Mass fraction n-heptane 
s = Product mass flow rate 


(+b). 
Fig 12. Comparison of centre-feed" operation 
with two-feed cocurrent operation. 
Legend : 
This investigation : 
TF-SCO-K, AT = 318 C 
TR-SCO-K, AT = 395°C 
a em 
Data of Lonemine [2] ; 
Centre-feed operation 
+ AT = 372°C; a=0079 cm 
OAT = 360°C; a = 0043 cm 
w = Mass fraction n-heptane in product 
P = Product mass flow rate. 


present study for two-feed cocurrent operation ; 
the separations were performed with a feed 
mixture of 0-500 mole fraction n-heptane in 
benzene at a A T for the upper curve of 39.5 C 
and for the lower of 31-8 C. The experimental 
points indicated for “ centre-fed ” operation are 
from Lonemire [2]; the separations were per- 
formed with an average feed concentration of 
0.500 mole fraction n-heptane in benzene. The 
upper of the two curves was obtained at a AT 
of 37-2 C and a plate spacing, a, of 0.079 em ; 
the lower of the two curves was obtained at a 
AT of 36-0 °C and a plate spacing of 0.043 em. 
The remaining column specifications for the two 
curves are the same; the thermogravitational 
column used by Loncmine was a concentric tube 
device with a height, h, of 183m and a heat 
transfer surface area of 961-8cm*. By way of 
comparison the experimental data for two-feed 
cocurrent operation (this investigation) were 
obtained with a flat plate thermogravitational 
column with rservoirs; the column had a heat 
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transfer area of 779-4cem? with a height, h, of 
19-33 cm and a plate spacing of 0.0376 cm. 

For the particular conditions described above it 
can be seen that two-feed cocurrent operation as 
carried out in this investigation produces much 
larger separations at the higher flow rates than 
obtainable from centre-fed’ operation. The 
column with a plate spacing of 0-079 em used in 
obtaining the upper of the two “ centre-fed ~ 
curves gave a batch steady-state separation of 
0-376; the other centre-fed” column with a 
plate spacing of 0.043cem gave a batch steady- 
state separation of 0-961. The column of this 
investigation gave a batch steady-state separation 
of 0-477. These separations are all mass fraction 
differences. 


Two-feed countercurrent operation 


Figure 13 shows the experimental results for 
“symmetrical” two-feed countercurrent opera- 
tion with a feed mixture composition of 0-150 


100 
(1/5), 
Fic. 13. Upper and lower reservoir product concen- 
tration difference vs. reciprocal average product rate : 
two-feed countercurrent operation, 
Legend : 
Two-feed countercurrent experimental point : 
TF-SCC-D, 1-6; AT = 821°C 
wp 0150 
- Two-feed cocurrent experimental curve 
Two-feed countercurrent equation (20), 
a 
———— Two-feed cocurrent, equation (15) 
w = Mass fraction n-heptane in product 
Arithmetic average product mass flow rate 
= + 


mass fraction n-heptane; the column was operated 
at a AT of 32:1 °C. The experimental points 
agree quite well with the predicted curve calcu- 
lated from equation (20) at the higher flow rates ; 
the two curves are not in agreement for the larger 
separations as a result of assumptions involved 
in the derivation of the separation equation. 

Also shown in Fig. 13 are the experimental and 
theoretical curves for “symmetric two-feed 
cocurrent operation taken at approximately the 
same operating conditions as the data for two-feed 
countereurrent operation. The cocurrent and 
countercurrent curves coincide at the higher 
product mass flow rates; at lower product flow 
rates two-feed countercurrent operation gives 
larger separations than found for cocurrent opera- 
tion. The theoretical curves for the two methods 
are also shown and substantiate the experimental 
results. 

Figure 14 presents experimental results for two- 
feed countercurrent operation with feed mixture 
composition of 0.500 mole fraction n-heptane 
in benzene and at a AT of 32.5 C. 


i 1 
100 
(i/5), sec/g 


Fic. 14. Upper and lower reservoir concentration 
difference vs. reciprocal average product rate: two- 
feed countercurrent operation. 

Legend : 
© Two-feed countercurrent experimental point : 

runs TF-SCC-K, 1-4 
() Two-feed cocurrent experimental point: runs 
TF-SCO-K, 1-4, 6 
a2 = Mole fraction n-heptane in product 
3 = Arithmetic average product mass flow rate 
=4(s, + 
Chem. Engng. Sci. Vol. 13, No. 4 June, 1961. 
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Single-feed cocurrent operation 

A difficulty with the theoretical discussion of 
single-feed cocurrent and countercurrent operation 
is that the modification of the “ transport ” 
equation (22), introduced the restriction that the 
flow parameter, r, be small so that the free and 
forced convection flows may be superimposed. 
This restriction is in conflict with the approxima- 
tions of equations (27) and (28); therefore, the 
area of agreement between experimental and pre- 
dicted results for single-feed operation should be 
smaller than found with two-feed operation. 

Figure 15 presents experimental separation data 
for single-feed cocurrent operation at a feed 
concentration of 0-500 mole fraction n-heptane 
in benzene; the operating AT was 32-4 C. 
Other operating conditions for the data were the 
same as described in previous experimental work. 


Fic. 15. Upper and lower reservoir product concen- 
tration difference vs. sum reciprocal product rates : 
single-feed cceurrent operation. 

Legend : 

© Single-feed cocurrent operation : 
Runs SF SCO_K, 1-4 
Two-feed cocurrent experimental curve, K 
Single-feed cocurrent, equation (39) 
a = Mole fraction n-heptane in product 
s = Product mass flow rate 


The theoretical curve calculated from equation 
(29) predicts the shape of the experimental curve 
but does not provide a quantitative prediction. 
A comparison between two-feed and single-feed 
cocurrent operation is shown in the Fig.; the 


two-feed cocurrent experimental curve taken at 
approximately the same condition is plotted, 
Experimental data for single-feed cocurrent 
operation and theoretical separations, calculated 
from equation (30), are given in Fig. 16 for a feed 
concentration of 0-150 mass fraction n-heptane 
in benzene and a column A T of 32-1 C. Also 
shown in this Fig. is the experimental two-feed 
cocurrent curve taken at the same conditions. 


00 
(1/5), sec’ 


Fic. 16, Upper and lower reservoir product con- 
centration difference vs. reciprocal average product 
rate ; single-feed cocurrent operation. 
Legend : 
© Single-feed cocurrent experimental point: Runs 

SF-SCO-D, 1-5; AT = 82:1 °C 
~~ Two-feed cocurrent experimental curve, 
TF SCO_-K 
~~ Single-feed cocurrent equation (30) 
w = Mass fraction n-heptane in product 
§ = Arithmetic average product mass flow rate 
(s, 


Single-feed countercurrent operation 


Figure 17 shows the experimental separations for 
“symmetric” single-feed countercurrent operation; 
the experimental data were taken at a feed 
concentration of 0-150 mass fraction n-heptane in 
benzene and a column AT of 32-4°C. Other 
conditions are the same as described for previous 
experimental work, Single-feed countercurrent 
separations are compared with two-feed cocurrent 
separations taken at the same conditions ; 
two-feed cocurrent operation gives greater separa- 
tions at the higher flow rates, but at the lower 
flow rates single-feed countercurrent separations 
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(45), sec’ 


Fic. 17. Upper and lower reservoir product con- 
centration difference vs. reciprocal average product 
rate: single-feed countercurrent operation. 
Legend : 
© Single-feed countercurrent experimental point : 

Runs SF-SCC-D, 1-6, 8; AT 32-4 °C 
~~~ Two-feed cocurrent experimental curve, 
TF-SCO-D 
Single-feed countercurrent ,equation (30) 
w — Mass fraction n-heptane in product 


are greater, Just as for two-feed operation, the 
single-feed cocurrent and countercurrent experi- 
mental curves coincide (this coincidence is not 
shown) at the higher product flow rates. At lower 
product flow rates the experimental separations 
for single-feed cocurrent operation become greater 
than those for single-feed countercurrent opera- 
tion ; the opposite result was found for two-feed 
operation. 

Also given in Fig. 17 is a comparison between 
the experimental separations found for single-feed 
countercurrent operation and the separations 
predicted by equation (30), The agreement is far 
from quantitative for reasons mentioned earlier. 
However, the theory correctly predicts the 
experimentally observed coincidence of the sep- 
aration curves for single-feed cocurrent and 
countercurrent operation at the higher flow rates. 

Experimental data for “ asymmetric” single- 
feed countercurrent operation are given in Fig, 18 
for feed mixtures of 0-150 mass fraction and 0.500 
mole fraction n-heptane in benzene. For these 
runs the upper reservoir product flow rate is fixed 
and the lower reservoir product flow rate is varied. 
The data for the dilute case indicate that the lower 


. 


*(0-01323,0034) 


*(0-0221,0:020) 


0-002 0-004 0-006 


S, g/sec 


Fic. 18. Feed and lower reservoir product con- 
centration difference vs. lower reservoir product flow 
rate: single-feed countercurrent operation. 


Legend : 

© & = 00231 g/ sec 
AT = 820°C 

Runs SF-USCC_D, 6-8, 3 
0-02606 g sec 
= $2-1 °C 

Runs SF-USCC-K, 1-3 
s, = 0-0231 g ‘see 
AT VC 

Runs SF-USCO-D, 1, 2 

w = Mass fraction n-heptane in product 
s = Product mass flow rate 


AT 


reservoir product becomes pure benzene at a 
finite flow rate; this result suggests that this 
method of operation could be quite useful for 
the purification of substances. 


CONCLUSIONS 


1. Theory and experiment for two-feed and 
single-feed operation are in qualitative agreement ; 
quantitative agreement is obtained for the two- 
feed operation methods at the higher product 
flow rates through the use of experimentally 
determined “ transport constants.” 

2. Two-feed cocurrent and countercurrent 
operation of a flat plate thermogravitational 
column with reservoirs should provide additional 
methods for the experimental measurement of the 
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liquid phase “ transport constants’ and refine- 
ment of thermogravitational column theory. 

3. All four of the operation methods discussed 
may be of use in liquid mixture separation 


problems. 
4. It is suggested that the two-feed operation 


methods coupled with a method for controlling 
the natural convection (applied electric fields or 
moving walls would be two possible methods) 
would provide a relatively eflicient method for 
effecting separations in the liquid phase with 
thermal diffusion. 
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NOTATION 


Perpendicular distance between “hot” and 
“cold” walls in a flat plate thermogravitational 
column 
Slope of linear approximation to @,@, for dilute 
binary mixtures 
Intercept of linear approximation to @,@, for 
dilute binary mixtures 
Mutual diffusion coefficient in binary mixtures 
Transport constant from solution of continuity 
equations for the physical system of this investiga- 
tion 
Vertical height of thermogravitational column 
Mass flux vector with respect to mass average 
velocity 
“ Transport constant " from solution of continuity 
equations for the physical system of this investiga- 
tion 
Horizontal length of flat plate thermogravitational 
column 
Net rate at which component “ 1" is moving up 
the column 
Net rate at which component “1” 
down the column 
Product mass flow rate 


is moving 


Feed mass flow rate in single-feed operation 
Arithmetic mean temperature of “ hot" and 
cold” walls 
Mass average velocity 
Cartesian co-ordinate perpendicular to heat trans- 
fer surface 
Mole fraction of i-th component 
Cartesian co-ordinate parallel to heat transfer 
surface 
Cartesian co-ordinate parallel to heat transfer 
surface, horizontal direction 
Volumetric coefficient of thermal expansion 
Angle of flat plates from the vertical 
Coefficient of shear viscosity 
Mass density 
Soret coeflicient of component “1” 
Flux of * 1” in z direction produced by concentra- 
tion diffusion 
In a binary system the mass fraction of that 
component which in the presence of a temperature 
gradient migrates to the “ hot” wall and con- 
centrates at the top of a thermal diffusion column 
@, Mass fraction averaged over x, dz horizontal plane 
at a vertical position, y 
F Quantity pertaining to feed 
(subseript) 
i =Quantity pertaining to lower reservoir in thermo- 
(subseript) gravitational column 
s Quantity pertaining to upper reservoir in thermo- 
(subscript) gravitational column 
Vector 
(above symbol) 
Explanation of code in Figs. 
K Feed composition 6-500 mole fraction n-heptane 
in benzene 
D Feed composition 0-150 mass fraction n-heptane 
in benzene 
S “Symmetrical operation ; equal product mass 
flow rates 
Asymmetrical" operation; unequal product 
mass flow rates 
Two-feed operation 
Single-feed operation 
CO Cocurrent operation 
- Countercurrent operation 
Column data for all runs 
a = 0-0376 em 
h = 19-33 em 
L = 40-32 on 
O = G1 deg. 


Us 


TF 
SF 
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Intraparticle diffusion and conduction in porous catalysts—I 
Single reactions 


R. E. Scumson* and Neat R. Awunpson 
University of Minnesota, Minneapolis 14, Minnesota, U.S.A. 


Abstract Calculations are made for a single porous catalytic sphere in which a single chemical 


reaction is carried out. Transport within the sphere is by Knudsen diffusion. Heat is carried by 
molecular flow and by conduction through the solid structure. Approximate and rigorous caleula- 
tions are made for partial pressure and temperature profiles and an approximate solution based 
on a linearization of the heat generation function is shown to be adequate for a number of cases. 
Effectiveness factors are computed based on the external surface conditions. 


Résumé— Cet article traite des calculs relatifs & une réaction chimique simple se produisant au 
contact dun catalyseur poreux de forme sphérique. Le transfert dans la sphére se fait par diffusion 


de Knudsen. 


L*écoulement moléculaire et la conduction a travers la structure solide assurent 


le transfert de chaleur. Des caleuls approximatifs et rigoureux de la pression partielle et des 
profils de température ct une solution approximative supposant la linéarité du dégagement, 
conviennent dans la plupart des cas. Les facteurs d’activité sont calcules d'aprés les conditions 


de surface extéricure. 


Zusammenfassung 


cinzige chemische Reaktion stattfindet, sind durchgefiihrt. 


Berechnungen fiir eine cinzige katalytische Kugel, in welcher sich eine 


Beforderung innerhalb der Kugel 


ist durch Knudsen Diffusion. Die Wirme ist durch molekularen Fluss und Leitung durch die solide 
Struktur tibertragen. Annihernde und genaue Berechnungen fiir partielle Druck- und Temper- 
atur-profile sind gemacht, und es ist gezeigt, dass cine annihernde Losung, auf der Grundlage der 
Linearisicrung der Wirmeerzeugung-funktion in zahlreichen Fallen hinreichend ist. Wirksam- 


Tue purpose of this paper is to discuss and develop 
methods for the prediction of the rate of reaction 
from a single porous catalyst particle under the 
assumption that the kinetics of the reaction are 
known, that Knudsen diffusion obtains in the 
void volume of the catalyst and that the only 
mode of heat release is by intraparticle conduction 
in the porous structure. Simpler problems have 
been considered and the main deterrent to progress 
has been the fact that the equations which need 
solving are highly non-linear. In fact, the only 
system which yields analytic solutions is the 
simple reversible reaction first order in each 
direction for a reaction carried out in particles of 
spherical, cylindrical, plate-like, or rectangular 
catalyst shapes at constant temperature. Prob- 
lems of this kind were considered by Tarene [1] 
and Ze.powrrscn [2] almost simultaneously. 
Wueecer [3] and Prarer and Wertsz [4] have 


keits-faktoren sind den dusseren Oberflichen-verhiltnissen ertsprechend berechnet. 


*Present address: The Ohio Oil Company, Littleton, Colorado, U.S.A. 


discussed the problem of intraparticle diffusion 
in some detail generally for isothermal operation. 
The former, along with Arts [5] has considered 
the effect of geometry of the catalyst while 
Prarer and Weisz have made use of Tue e’s 
and Wuerr.er’s results to interpret data and 
design experiments. Finally, Wicker [6-9] in a 
number of papers has discussed the effect of 
intraparticle diffusion in carbon burning. Srur 
[10] has investigated thermal conductivities of 
catalyst pellets. Many of the pertinent literature 
references in the field are given by Wetsz in 
Ref. [11] and others [3, 12-18}. 


EQUATIONS FOR A SINGLE REACTION 
n 
Consider a single reaction 2 a;M,; = 0, a; > 0 
int 


for a product and a; < 0 for a reactant, taking 
place in a porous catalyst pellet of spherical shape. 


Suppose the particle has randomly interconnected 
pores and a fractional void volume y. It is 
assumed that reactant and product molecules 
diffuse through the fractional void volume. The 
reaction rate within the sphere is given by f,; where 
J, is the rate of formation of component # in 
moles per unit time per unit area of catalyst 
surface. If a mass balance on a spherical shell of 
inside radius p and thickness Ap is made, then 
the rate of diffusion through the inner surface 
outward is 

(4x9? dD, dp, ‘) 

dp ° 

in moles per unit time where the subscript p 
indicates that the quantity in parentheses must 
be evaluated at p. D, is the diffusion coefficient 
computed from 


4 Ss E, 
2 By2r7RmT T 
in which the factor y/2 has been suggested by 


Wherever. The rate at which component / leaves 
the spherical surface at p + Ap is 


(tp 2D, 
‘dp 
and the rate at which M, is formed in the reaction 
in the volume 42 p* Ap is 
Ap (47 py 


where p < p < p + Ap. S, is the surface area per 
gram of pellet and p, is the pellet density. Sum- 
ming these three quantities, dividing by Ap and 
passing to the limit, there results 


+ p, = 0, i= 1,2..." (1) 
From the stoichiometry of the reaction it is clear 


that 
Sj 


i 


so that for any other specie M; 
d ay 
dp (> D, iz) +P 
Appropriate manipulation of these two equations 
gives 
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1d Ld dp, 
a, dp \P a, dp \? “ 


Two successive integrations give 


E, E, 
— Pe) = (Pj Pjs) 
or 


P; = = + Pis) 


(=) (ri — (3) 


Thus instead of there -“ a differential equation 
for each component as in equation (1), there is 
a single differential equation for the conservation 
of mass, since in f,, which may be a function of 
the partial pressures of each component, one 
may use equation (3) to eliminate all partial 
pressures save that of the i-th. f, is then a function 
of one partial pressure p, and all partial pressures 
at the surface of the particle, i.e. 


Si = (Pir Pie Paw Pus T) 

Let us consider now the equation of conservation 
of energy within the particle. It will be assumed 
that the diffusing molecules carry their partial 
molar enthalpy, and, since only Knudsen flow 
is being considered, this is the only mode of 
transfer in the gas phase. In the solid phase heat 
will be transported by conduction. It should be 
noted that the partial molar enthalpy should be 
corrected for the heat of transport. Once again, 
analysing a spherical shell, the rate at which 
heat enters the inner surface is 


te ((1 — y)k, p? = 
dp 
by intraparticle conduction and 


by molecular flow. 
It may then be shown that the conservation 


equation is 
+ — A,| =0 
One integration may be easily made 
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_D; dp, 
4; dp 


This when substituted into equation (4) gives 


a; dp 
Now it may be shown easily that in molecular 
flow through a capillary A; is not the partial 
molar enthalpy but rather 71, the partial molar 
enthalpy diminished by 4 RT, so 


— i, 
Thus 
= Fa, 


tel i=l 
AH, —4 RT Za, 
i=l 


where AH, is the heat of reaction of the reaction 
n 
z a; M, = 0. 
Therefore 


(1 — yk, 


La) =0 

It will be siaeail that the term involving 4 RT 

is small compared to AH, so that 


a 


This equation may be used to estimate the 
maximum temperature difference within the 
particle if it is assumed that p, falls to zero 
somewhere. 

Then 

pi, D, AH, 


(AT 
(l — y)k, aj 


Equation (5) may be rearranged and integrated 
to give 


= Py = TaT (6) 


and since k, and AH; may be expressed as funce- 


tions of the temperature the integral may be 
evaluated, Equation (5) may also be written 


(1 y) ky a; 
T Tf, 7 
where the parameters must take on average 
values over the temperature range T — T,. If k, 
and AH, are assumed to be constant equation (7) 
may be written 


Pj = Pis 


Pj = Pr — ( T°" — 
and this along with equation (7) will be used in 
the examples. 
From equations (1) and (3) one may now write 


d (1 — y)k,a,dT 
dp dp 


where /,, as written in this formula, is a function 
of the temperature alone, the ith partial pressure 
having been eliminated through equation (6) or 
its modification. Equation (8) then is a fairly 
rigorous equation and its solution along with the 
boundary conditions 


aT 
dp 


Sy Si =O (8) 


0 at 


should give the temperature profile within a 
single catalyst particle. 

The assumption will now be made that hk, and 
AH, are constant so that equation (8) may be 
written in the form 


2dT 
dp?” p dp 
and with rR = p 


Pp” S, AH, 


9 


Rp, Sy 


0 (10 
(1 y) ky 


(11) 
(12) 
In addition to the temperature profile one is 


interested in the effectiveness factor for the 
catalyst. The effectiveness factor is defined as 
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dp int dp 
and from a single integration of equation (2) 
D, dp, 
a; dp 
Dz, dp; 
eT  2dT 
|| 
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the rate of production of a given chemical specie 
by a single particle divided by the fictitious rate 
assuming all the catalytic surface is available to 
reactants at the temperature and partial pressures 
at the surface. The rate of production per particle 


— dp, 


surface 
The fictitious rate is 


and the effectiveness factor is 


In terms of the eiainae profile 


pp R AH dr surface 


SOLUTION FOR A SINGLE REACTION 


The purpose of the remainder of the paper is 
to show how solutions of equations (10)-(12) may 
be obtained approximately and rigorously, the 
approximate method by an analytical formula 
and the rigorous method by a technique which 
will give answers to any preassigned accuracy. 


Before proceeding, however, the function 
Pp . (— ¢ 
(1 — a, 


will be examined, Note first that ¢(7') has the 
sign of (— AH,) since f,a,;"' is always positive. 


q(T) = — 


In general it has been found that for a large 
number of cases the function q(7') takes on the 
general forms shown in Figs. I(a), 1(b) and 1(e). 
Figs. I(a) and 1(b) are typical of exothermic 
reactions while (Fig. 1(c) is typical of endothermic 
reactions. The temperature 7, is defined by 


q(T ,) 


and gives the temperature at which either the 
limiting reactant has been exhausted or is, alter- 
natively, the point at which chemical equilibrium 
has been established in the particle. In any 
event T, — T, is a measure of the maximum 
temperature difference in the particle and so is a 
generalization of (AT) previously mentioned. 

First an approximate solution of equation (10) 
will be developed, Since it is apparent that the 
heat generation function g (7) may be approxim- 
ated in two of the cases, Figs. I(a) and Fig.1(c), 
by a straight line with good accuracy, the temper- 
ature profile with a — bT replacing q(T) should 
give a good approximation to the real temperature 
profile, Thus 


2 dT 


te 


needs solving. This is a straightforward problem 
and the solution may be written 


T (a, sinh r b 


T (r) - ; +o! (13) 


Sinh yb 


(a) 


(do) (c) 


Heat generation function, q(T) 


Temperature, T 


Temperoture, T 


Temperature, T 


Fic. 1. Typical heat generation function. 
(a) Exothermic reaction type 1 (b) Exothermic reaction type 2 (c) Endothermic reaction 
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To solve equation (10) the following procedure 
will be adopted. Write it in the form 


_2 dT 


dr? sor dr 
then 


+q(T)+a—bT =a—bT 


@T 2adT 

where 


Q(r) a —bT —q(T) 


If one supposes for the moment that Q(r) —a 
is a known function of r then the solution of 
equation (15) may be generated by methods 
readily available to give 

T, — (a/b) sinhr yb 
sinh 

1 
+ [ ee. (t) dt (15) 


T(r) 


bsinh 


O<xt<cr 
sinh — t)sinh y br 


ry bsinh yb 


r<xt<l 
The solution, equation (15) consists of three terms. 
The first two terms on the right-hand side make 
up the approximate solution. The third or integral 
term may be considered to be a correction term 
to this first approximation. In order to evaluate 
the correction term a method of successive 
substitutions will be applied. Q(t) depends upon 
T (r), the temperature profile, which is at this 
point unknown and therefore equation (15) is 
an integral equation with the unknown function 
T (r) appearing on the left as well as under the 
integral sign. If the function q(T) is almost 
linear q(T) — a +-bT is small provided that a 
and b are chosen properly. If the process of suc- 
cessive substitutions converges, then using the 
temperature profile, T(r), the profile obtained by 
assuming q (7) = a — bT asa first approximation 


Q, (t) = a — 6T, —q(T)) 


then 
1 
T,(r) = T,(r) + 


should be a better approximation. The process 
may be repeated using Q,(r) — a — bT, — q(T,) 
to define a new profile 

1 


T;(r) = T,(r) 4 [ 


and in general 
1 


with ,(t) — a — bT,_, —q(T,,,), and so a 
sequence of profiles T,, T,, Ts... . ete. may be 
computed to give ultimately a profile 7',(r) which 
differs from the previous profile 7),,(r) or the 
succeeding profile 7,,.,(r) to any preassigned 
accuracy. For Type Ma) and Type I(c) heat 
generation functions this process appears to 
converge and produces profiles which satisfy the 
differential equation and boundary conditions 
numerically. 

A word is needed concerning the method of 
evaluation of the integrals. Many methods for 
the numerical evaluation of integrals may be 
suitable but in this paper the integral was replaced 
by a sum, using the trapezoidal rule for each 
sub-interval. The interval from zero to one was 
broken into twenty equal subintervals. Greater 
accuracy may be obtained with a finer subdivision 
but in the cases investigated here, the number 
twenty appeared adequate. The choice of which 
linear function a — bT to be used is not difficult, 
One may use the Taylor series approximation 
about the surface conditions or one might fair 
a straight line to q(7') by eye as shown in Fig. 2. 

If the heat generation function q(T’) is of Type 
I(b) the method described above may be un- 
satisfactory, but as shown in Fig. (5) the curve 
may be approximated by two straight lines. Let 
the point of intersection of these lines be r, and 
the temperature at r, be T,. The general problem 
to be solved then reduces to the set of equations 


VOL, 

with 7 

and 
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(18) 
(19) 

eT, 247 
t by Ty — ay (20) 
T,-T, (21) 
T,=T, (22) 

dT, dT, 

OL. dr dr 


with 
QA(r) = a, —b,T —q(T), 
= 49 + T — Q(T), 


The solution of this problem is obtained by 
solving for T,, from equations (17), (18) and (19) 
and for T, from equations (20), (21) and (22). 
These two solutions are obtained by a straight- 
forward method of variation of parameters and 
will involve r, and T', as 


O<r<r, 


, (a, ‘b,)] sinh r 


r sinh y/(b,) r, 


+ | G,(r.)tQ. (24) 
0 


_ sinh ¥/(6,) (r, — r) sinh y/(b,) ¢ 
r V(b.) sinh y/(b,) re 
Ostar 
sinh v r sinh (rz — 4) 
rv/(6,) sinh re 
r<t<r, 


sin (bo)(1—r) 
by! sin (bo) (1 — 


+ { G, (r, t) t Qy (t) dt (25) 


sin \ (1 —r) r)sin y (bo) (t 


ry (bo) sin y (by) — r,) 


r<t<r 
G, (r,t) =< 
sin y (bo) (r — r,) sin y (b9)(1 


(bo) sin y (bo) (1 — 

a r<t<l 
The temperature 7, may be determined from the 
solution of 


—6b,T, 


as a first 

In order to determine the value of r, the two 
equations above must be differentiated and evalu- 
ated at r= r,, 


=a +b, T, 


A first approximation to the temperature profile 
may be obtained from the equations 


Ta = 


T, — (a,/b,) sinh r y 6, 


O<r< 
r sinh r, yb, 


0 

4 r,(T, t | sin (b,) (1 —r) 
r sin (1 — r,) 


(27) 


where r, is determined by equating the derivatives 
atr=r,, or 


(7, coth yb. r, — — 


by 
T, do \ 


Further approximations may be obtained by the 
successive substitutions 


+ 


(28) 


To = Ta + 
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and 
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+7 ("+ = 
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with T,, and T,,, having obvious representations. 

It is apparent that a correction must be made 
to r, at some point. This need not be done at 
each iteration since it is found that r, does not 
change appreciably. It should be calculated, 
however, and it may be obtained by equating the 
derivatives of equations (24) and (25) and solving 
for r,. 


NUMERICAL EXAMPLES 


Four numerical examples will be given which 
should serve to illustrate the principles developed 
above. The diffusion coefficients will be calculated 
from 

59-109 

Example I. 24 +B 
AH, = — 140,000 cal/g mole B R — 0.4762 em 

m, 45, mp — 90 385A 

T, 600 °K k, = 8 10-4 

Pa, 3-20 atm cal ‘sec om °K 

y= 0-35 
ta 2-888 (10 *) exp 9058 T) p,?. 
g mole A sec 

D, = 4407 g mole atm see em 
q(T) = 4-703 (10-*) p Pexp(— 9058 T), “°K 

P4 = 41-13 — 0.06321 T, atm 


Fic. 2. Heat generation function for Example I 
and two linear approximations. 
1- — a — bT = 15487 — 2458T 
2 bT = 18480 — 2-57T 


Fig. 2 gives the graph of ¢(7') vs. T and shows 
two linear approximations, 


line 1 ~a —bT — 15,480 — 23.58T 
line 2~a — bT = 18,480 — 28.577 °K 


Fig. 3 shows the temperature profile within the 
particle using line 1 as a first approximation. 
The second and third approximations are also 
shown using equation (16) as well as the higher 
order approximation assumed to be the rigorous 
solution. Fig. 4 shows the first approximations 
using lines 1 and 2. The effectiveness factors are 


0-461 
0-490 
0-456 


rigorous solution 
line 1 approximation 
line 2 approximation 


Example Il, A +B 
AH, = — 30,000 cal/gmole = 0-5em 
m, = 36 y 
= 625 °K y 0-40 
atm k, 80 10-4 
cal em °K 
0.06254 exp (15,000 p mole 
A em? sec 
115-74 — 6-767 (10°) T? 2, atm 
== 19542 (10") exp (— 15,000/T) py. “K 


660 


| 
02 06 
Dimensioniess rodius, 


Fic. 3. Temperature profile for Example 1. 
first: approximation 


second approximation 
— third approximation 
~ final approximation 


= 
§ 620 
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02 O4 Ot 


13 Dimensionless radius, 


Fic. 4. Comparison of true and approximate temp- 
erature profiles for Example 1, 

Final solution, rigorous method 

~ First approximation, line 1 
First approximation, line 2 


Fig. 5 shows the graph of g (7') vs. the temperature 
and two separate approximations, one by a single 
line and another by two straight lines 
one line ~a — bT = 12,620 — 19-01T, K 
two lines ~ a, — b, T = 25,530 — 38-467, °K 

dy + by T = — 2194 + 4-691T °K 
Fig. 6 gives the temperature profiles calculated 
from the rigorous solution and the two first 


62 


Temperoture, T, 
Fic. 5. Heat generation functions for Example I 
showing one- and two-line approximations. 
« bT Linear approximations, two-line method 
a—bT Linear approximation, one-line method 


665 


Temperature, Tir), 


Dimensionless radius, 


Fic. 6. Temperature profiles for Example II. 
Final solution 

First approximation, two-line method 

— First approximation, one-line method 


approximations, one using a one line approxima- 
tion for g (7) and the other for a two-line approx- 
mation. It should be pointed out that even 
though the single straight line appears to be a 
poor approximation the final results by the two 
methods are essentially identical although twice 
as many successive substitutions are required in 
one than in the other. 


One line 
First approximation E.F. 0712 0-531 
Rigorous solution E.F. 0-709 0-713 
Centre temperature 663-6 663-6 


Two lines 


In the two-line method the first approximation 
is almost as good as the third approximation, 
while the third differs almost imperceptibly from 
the ninth. 


Example 24 +B 
AH, = — 160,000 cal /g mole B R = 0-25 em 
m, = 36 ¥=200A 
T, = 400 °K y = 0-400 
Pa, 10 atm k, = 8 x lo-* 
cal /sec em °K 
0-702 (10"") exp (— 20,000 /T) p_,?, 
g mole A /cm? sec 
3-941 (10-7), g mole /atm cm sec (assumed 
constant) 
P4 = 7-090 — 0.01523 T , atm 
q(T) = 2-925 (104) exp (— 20,000/T') p 4?, 
Chem. Engng. Sci. Vol. 13, No. 4. June, 1961. 
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“« 


generction function, 


a0 
Temperature, T, 


Fic. 7. Heat generation function and two-line 
approximation for Example III. 


approximations, two-line method 


a, -&T 


Fig. 7 gives the graph of g (T') vs T, and it appears 
that q(T) can only be approximated by two 
straight lines 


a, 6 T 22,868 — 49.107, 


ao by T 6706.5 + 16-91 
and the intersection point is T, — 448 °K. The 
distance r, can be computed for the first approx- 


470 


Dimensioniess radius, 
Fic, 8, True and approximate temperature profiles 

in Example ITI. 
Final solution 
~---~- First approximation 


imation and is 0.665, After a few approximations 
r, was recalculated and was found to be 0-672. 
From the first approximation the effectiveness 
factor was 7-65 while for the rigorous solution 
it was 7-29. Fig. 8 shows how closely the tem- 
perature profiles of the first approximation and 


the rigorous solution agree. 
Example IV. 24 +B 
All; — 160,000 cal g mole B 

ma 

T, 600 °K 

Pa, = 2-0 atm lot 
cal sec cm “K 
10,000, T) p 


ta 0.775 (10 exp ( 
g mole A em? sec 
1.931 + 2-675 (10 4) 7°? atm 


2-583 (10%) exp (— 10,000/T) 


Pa 
q(T) 


Temperature, T, °K 


Fic. 9. Heat generation function and two approx- 
imations for Example IV. 


Linear approximation, one-line method 
approximations, two-line method 


a—bT 
a, 
a,—b,T 


Fig. 9 gives q(T) vs. T and shows possible one 

and two line approximations. Suitable equations 

are 

one line = a — 86,180 —7025T, 

two lines ~ a, — b, T = 8209 —16-42T, 
dy — by T — 82,600 — 147-6T, °K 


O00 
V4 
j 
400 
VOL. 
13 
| | 
\ 
-2000}— 
> 500 $20 $40 S560 580 600 
| | | | 
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Dimensioniess radius, © 
Fre. 10. Temperature profiles in Example TV. 
Final solution 
First approximation, one-line method 
First approximation, two line method 


with r, — 0.889 and T, — 567 °K. Fig. 10 gives 
the usual comparison of the temperature profiles. 
The effectiveness factors are 0-241, 0-315 and 
0-240 for the rigorous solution, the one line first 
approximation, and the two line first approxima- 
tion, respectively. 


SUMMARY AND CRITIQUE 

Techniques have been presented which enable 
one to treat the problem of non-isothermal 
reaction within a porous catalyst pellet. The 
most important result is that the two approxima- 
tion methods presented, either the one-line or 
two-line methods, give results which indicate 
that for single reactions, at least, a more detailed 
solution is not essential. The pertinent point here 
is that the heat generation function, ¢ (7'), must 
be examined and it was apparent that in spite 
of the fact that it appeared to be rather complex 
it could here be satisfactorily approximated to 
in the four numerical examples presented. 
Problems with more complicated kinetic expres- 
sions involving adsorption denominators have been 
investigated and as yet no case has been con- 
fronted which could not be treated by one of 
the two techniques. 

The problem with a film resistance at the 
boundary may be treated in the same way 
except that the details of the computations differ. 
Obviously, the problem of the isothermal particle 


with complicated chemical kinetics follows directly 
also from the above analysis. 

One point which needs some elucidation is the 
method of numerical solution for what has been 
called above the rigorous solution, although one 
result of this paper is that such calculations 
need not be made. In order to solve the different- 
tial equation (10), it was changed to an integral 
equation but first subtracting out a first approx- 
ation. The integral in the integral equation was 
then changed to a sum and the set of nonlinear 
transcendental equations was solved by successive 
substitutions. Although no formal proof was 
presented the successive substitutions seemed to 
converge and it was assumed that the desired 
temperature protile was obtained. Now there 
are many other methods of solving such a set of 
equations ; however, many that were tried were 
not successful and some diverged badly even with 
a good first approximation. 

When one is presented with an ordinary 
differential equation today the analogue computer 
immediately comes to mind and this was tried 
here, albeit not exhaustively. These problems 
are two-point boundary value problems and so 
in order to use the analogue computer the value 
of the derivative must be estimated at the particle 
surface or the temperature must be estimated at 
the centre of the particle. While conceptually 
this is very simple, there are difficulties. It is 
found that the equations are so unstable numer- 
ically that there is no value of the derivative at 
the particle surface which will enable the equations 
to be integrated by a marching technique. Solu- 
tion on the analogue computer starting at the 
centre of the particle where the gradient is zero 
may be accomplished in some cases. With this 
procedure the temperature must be estimated at 
the centre and revised until the solution agrees 
at the surface with the prescribed surface tem- 
perature. This phenomenon is not peculiar to 
these problems, but from the experience of the 
writers, it does occur in many problems which 
involve combined chemical reaction, heat release 
and diffusion. The mathematical problem is 
extensively discussed in the book by Fox [19]. 

One further point needs discussion. The defini- 
tion of the effectiveness factor given by Taree 
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implies that it will always lie between zero and 
one. In problem III of this paper an E.F. of 
over seven was obtained and it is apparent that 
this is possible since the high temperature within 
the particle produces an abnormally high rate of 
reaction. 
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NOTATION 


constants in linearization of q (7) 

constants in linearization of q(T), outer section 
constants in linearization of g (7), centre section 
stoichiometric coeflicients 

Knudsen diffusivity multiplied by y/2 
temperature independent part of D; 

rate of production of M; in moles per unit area 
per unit time 

Green's function in solution for temperature 
profile 

heat carried by M,; per mole 

partial molar enthalpy of M; 


AHy 
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Abstract—The problem of intraparticle diffusion and conduction is considered for complex 
reactions taking place in porous catalytic spheres. It is assumed that Knudsen diffusion prevails 
and that heat generated in the reaction is conducted through the porous structure to the catalyst 
surface. Equations are written for an arbitrary number of reactions and criteria are developed to 
determine the number of conservation equations for complex systems. A particular numerical 
example is worked both approximately and exactly by the same technique as used in the previous 
paper; the term exactly meaning that partial pressure and temperature profiles may be calculated 
numerically as accurately as desired. The approximate method should be satisfactory for almost 
all industrial work. 


Résumé — Le probleme de la diffusion et de la conduction a Vintérieur d'une particule est considéré 
pour des réactions complexes se produisant avec des catalyseurs poreux de forme sphérique. On 
suppose que la diffusion de Knudsen prévaut et que la chaleur dégagée par la réaction est conduite 
jusqu’a la surface du catalyseur a travers la structure poreuse de celui-ci. Des équations sont 
posées pour un nombre arbitraire de réactions et des critéria développés pour déterminer le nombre 
d’équations de conservation pour des systémes complexes. Un exemple numérique particulier est 
étudié a la fois approximativement et exactement par la technique utilisée dans les précédents 
articles. Le terme exact montre que la pression partielle et le profil de température peuvent étre 
calculés numériquement de fagon aussi précise que lon veut. La méthode approximative devrait 
étre satisfaisante pour la plupart des cas industriels. 


Zusammenfassung-— Das Problem der intrapartikuliiren Diffusion und Leitung wird fiir 
komplexe Reaktoinen betrachtet, die in pordsen Katalysatorkugeln stattfinden. Dabei wird 
angenommen, dass Knudsen-Diffusion vorherrscht und dass die Reaktionswirme durch die 
porése Struktur zur Katalysatoroberfliche geleitet wird. Fiir eine willkiirliche Zahl von Reak- 
tionen werden Gleichungen angegeben und Kriterien entwickelt, um die Zahl der Bilanzgleichungen 
fir komplexe Systeme zu bestimmen. Ein numerisches Beispicl ist sowohl angenihert als auch 
exakt mit der gleichen Methode durchgerechnet, die in der vorhergehenden Arbeit benutzt wurde ; 
mit Hilfe des exakten Ausdrucks kénnen Partialdruck und Temperaturprofil numerisch beliebig 
genau berechnet werden. Die Niherungsmethode sollte fiir die meisten industriellen Arbeiten 


geniigen. 


arbitrary system of chemical reactions will be 
considered [2]. 

The chemical species will be denoted by M, 
and the stoichiometric coeflicient of the ith species 
in the jth reaction will be a;;. Suppose there are 


INTRODUCTION 


IN A previous paper [1] the authors considered 
the problem of a chemical reaction taking place 
inside a single porous catalyst sphere. It was 
shown how the temperature and partial pressure : 
profile could be calculated and how the effective. 1” Teactions 


ness factor could be estimated. An approximate 
method was developed and it was shown that the 
approximation made by a suitable linearization 
gave results of satisfactory accuracy. It is the 
purpose of this paper to generalize the methods 
of the previous paper to complex reactions. An 


2 ay M; = 0; j=1,3,...,@ 

i=! 
where in a particular reaction a, is positive for 
a product and negative for a reactant. The heat 
of reaction for the jth reaction is AH; and f;; is 
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the rate of formation of the ith species in the jth 
reaction in moles per unit of catalytic surface 
area per unit of time. The quantity f, may be 
a function of all x partial pressures, the temper- 
ature and the total pressure. We suppose that 
the set of chemical species {M;} contains all of 
the species in the reaction system. 

It is proposed to show how the partial pressure 
and temperature profile may be computed both 
exactly and approximately. The same methods 
may be used if the conductivity of the particles 
is great enough so that temperature gradients 
may be neglected. 


DERIVATION OF EQUATIONS* 


It is assumed that Knudsen diffusion obtains 
within the particles so that the diffusion coefficient 
is given by 


for the ith chemical species. E; is independent 
of the temperature and depends on the pore 
properties as well as the molecular weight. 
Assuming that all flow is by molecular diffusion 
a mass balance on a single particle is given by 


dj. 
D, =} — — 
z( ‘dp Pr Se 
228 


From the stoichiometry of the reaction it follows 
that 


Sig Sug 


1 for all j 
ay kj 


so that 
d ( +.) = 


— p* p, S, é=1,3,3,...9 (1) 


The question now arises as to the number of 
these equations which are independent. This set 
of equations may be written in vector form as 


dp (P29) —p*p, S,AT 
where q and f are vectors, 


- * Notation will be the same as in the first paper. 


and A is the matrix 


A=] ay 


Gace 


Consider the matrix A’, the transpose of A 


Ay, any 


AT 


and suppose A has rank s. For the sake of con- 
venience, suppose the nonzero determinant of 
order s is in the upper left hand corner. If it 
is not, the species may be renumbered to make 
it so. Now if the columns of AT are considered 
as vectors it follows that each of the last n-s 
columns may be expressed as a linear combination 
of the first s columns. Therefore, there is a matrix 
BY of the form 


O 0....0 Yn,2 


d 
dD, f 
d 
q dD, ; t= 
dp, 
; 13 
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and BT may be obtained from AT by operations 
on rows above, as by the Simplex Method of 
linear programming. In short, there is a matrix 
operator # such that 


AT — 


where #* is the matrix operator which operates 
on the columns of A, Note that the last n-» 
columns of BT represent the coeflicients in the 
expansion of the last n-s columns of AT in terms 
of the first s columns of A’, the first s columns 
being considered as a basis in a vector space. 
If one now writes 
Af 

it is aoe” that this is a column vector with 
elements ay and hence each of the last 


n-s deapletat in the vector is a linear combination 
of the first s elements as 


m 
Sj = 


4 


jel 

or from equation (1) 

i, ‘dp. dp 
i=s+lton 


This may be integrated once to give 


where use has been made of the fact that the 
gradients are zero at the centre of the particle. 
Since E; = D, /T the latter equation may be 


integrated to give 


(py 


— Pis) -2 Yiq E, (Py — Pa) 
i=s+lton (3) 


Therefore inside the particle there are s independent 
partial pressures since the remaining (n—s) may 
be eliminated from the rate expressions, and hence 
also there are s equations of the type given by 
equation (1). 


Equation (2), of course, may be written in the 
form 
q — Bq’ 
where 


The essential result of this development is 
that if the rank of the stoichiometric matrix is s 
there are s independent mass conservation 
equations, the remaining partial pressures having 
been eliminated by equation (3). 

Now the heat balance must be considered. 
The same assumptions are made here as in the 
previous paper so that the energy conservation 
equation is formally the same 


Ala — y)k, , DA; =0 (5) 


p 


Here again Ai, is the partial molar enthalpy and 
the heat of transport (— } RT) will be neglected. 
One integration can be made directly to give 


(1 — y)k, + 
dp 


Define now a vector H in which A, is the ith 
element, then 


dT 


+H q=0 (6) 


A vector AH of the heats of reaction may be de- 
fined where the ith element is the heat of reaction 
of the ith reaction. Hence 


ATH = AH 
Also the reaction system itself may be written 
AT™M 


where M is the vector of chemical species. The 
operator operating as 


will give another set of chemical reactions where 
the heat of reaction vector will be 


? ATH = an’ = BH (7) 
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where a typical element of AH’ is 


AH,’ = A, + t Yur 


r=ltos 


and also where the elements A//,’ are a linear 
combination of the original AH, From equations 
(4), (6) and (7) one may write 

dT 


(1 — +H" BQ —o 
dp 


The transpose may be taken to give 


a — yk” BH 
dp 


If AH’ is assumed to be a constant vector, then 


= (1 y) k, T3/*) 
E, AH; (py — Pis) 
i=1 


or if D, is taken a as constant 


s 
(1 — (T = — =D, AH,’ 
i= 
and hence the temperature or one of the indepen- 
dent p,;'s can be eliminated from the rate expression. 
Hence of the original +1 conservation 
equations there are only s conservation equations 
which need integrating. Which equations in a 
particular case will be integrated depends on 
what is desired. In the example to follow in 
which AH; and D; are assumed temperature 
independent, a partial pressure and the temper- 
ature will be used as independent variables. 


EXaMPLe 
We will consider a reaction system 


A+BoC 
2. 2-D+E 


in which the reactions are irreversible and second 
order in the forward direction. The reaction is 
carried out in a porous pellet with the following 
information available. 


AH, = — 80,000 cal /g mole 
AH, = 60,000 cal mole 
25 ¢/g mole 
mp 36 mole 
me = 61 g/g mole 
Mp — 72g mole 
mp = 50g mole 
T, = 550°K 
Pay = 1-50 atm 
Pa, 1-50 atm 
Po. = O10 atm 
Pos = Pre = 0-0 atm 
em 
200 
0.40 
8 secem K 
4 B-892 (10-7) g mole 
- 3-244 (10-7) g mole 
= 2-492 (10°77) g mole 
2-293 (10-7) g mole 
Dy = 2-752 (10-7) g mole atm em see 


7500 le A 
Say 3:36 (10°) exp ( PsP, 


atm cm sec 
atm cm sec 
atm em sec 
atm em sec 


cm? sec 


2» g mole C 
pe 


Seg 1-480 (104) exp (- 


cm? sec 


In this problem the rank of the stiochiometric 
matrix is two so that there are two conservation 
equations. The detailed manipulations will not 
be given here but on the assumption that AH, &,, 
and the diffusivities are constant the following 
relations may be derived. 


2Dp 
p Ds 


Po dD, Pas + 


D,Au, 
Dp D, All, .) 
+ D4 AM, Pa 7) 

Dy AH," * Dy AH, 


(Pn. 
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The two conservation equations may be written 


d,p 2dp, . 
+w(T, py) = 90 
4 2dT 


+ q(T. py) 0 


2 


(T. Pa) 


The solution of this set of equations follows in 
principle the solution for a single equation except 
that it is more complicated. Our purpose here 
is to sketch the solution not to present the 
detailed technique. A particular example will be 
solved rather than a general one but it will be 
obvious what changes need be made for other 
cases, It will be recalled that a certain linear 
approximation was made and that will be done 
here. Let 

and 


bog + + by T 


be linear approximations tog (7, p4) and w(T.p 4) 
respectively. Then the solution of the equations 


2dT 
a + dog + + 4, T = 0 


dr? 


with 


2 


bog + + T = 0 


should give an approximaticn to the temperature 
and partial pressure profiles, say the first approx- 
imation. This solution may be found by direct 


methods and it is 
T,(r) — X(r) ¥(r) 


Pai (r) = BX (r) « ¥ (r) 


where 
be —2a 
X(r = 
— x) 
D4, boo) — T, — oo) sinh wr (8) 
pw? — x) sinh p r 
, boo — Bagg 
a) 
Pas ~ 
9 
+ A (3 — x) sinha r (9) 
with 


[—u + — 4v)]'/? 


[—u — — 4v)]'/? 


= Ay + by 


7 ay — 


B 
+ — ayo 
With these functions the profiles may be plotted 
and this will be done. 
As a generalization of the previous paper the 
two differential equations are now written 


 2dT 
dr? + ra + + = Q(r) — 
4 


. 
bp, +byT = Wr) — 
r dr 


where 
Q (1) = + + — Pa) 
W (r) = bog + + On T — w(T. px) 


If it is assumed for the moment that Q(r) and 
W (r) are known functions of r, these two equations 
may be changed to integral equations by the 
method of variation of parameters to give 


T(r) X(r) ¥(r) + 


1 
+ | [xQ (e) — W («)] G,, (r, de — (10) 
0 
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Pal(r) = BX (r) — (r) + 


+ 8 W G, (r, e)de— (11) 


x | [AQ G, (r, ©) de 


Tesinh —r)sinhpe — 
— 
esinh pr sinh 


| (8 — 2) sinh 


and where G, (r, «) is the same except A is sub- 
stituted for ». Equations (10) and (11) are solved 
by successive substitutions to obtain a sequence 
of profiles which it is hoped will converge in the 
limit to the solution of the original set of differen- 
tial equations. That is, the profiles 7, (r) and 
Pa (r) are substituted into Q(«) and W(«) to 
obtain new profiles T,(r) and pyy(«) and this 
process is iterated to give a sequence of profiles 
{ T, (r), Pan (r)}. 

From these equations the effectiveness factor 
may be computed from 


KF —8D, (dp,/dr),., moles of M, 
S unit of time 
Pp 

A and » could be complex numbers but this alters 
the method in detail rather than principle. 

With the data presented at the beginning of 
this Section it may be shown that 

7500 

q(T. p4) = 5.00 (10%) exp (— 


Pa (11998 p, — 0-2997) — 


— 8-92 (10"*) exp 


(2-604 p, + — 39-12) 


w(T, p4) = — 8-682 (10%) exp 


Pa (1-1998 p, — 0-2997) 


where use has been made of the formula 
Py Sg = 2A/r. In order to study the nature of 
the functions g and w calculations were made and 
three-dimensional models made which are shown 
in Figs. 1 and 2. These models are interesting 
but they must be analysed further to find the 
region in the T.p, plane of interest. Since 
q(T,. p4,) is positive it appears that there will be 
a net exothermic effect so that T > 550° K in 
the particle. Since py, and p, must not be negative 


T > 609-2 — 40.55 


In addition py < 1-50. A further boundary line 
may be obtained from the condition that if the 
first reaction goes to completion and the second 
one is insignificant, then pp < 1-50 and therefore 


T 68246 — 0.55 py 


The resulting region of interest is shown in Fig. 
%. With this information the linear approxima- 
tions can be found and in this problem these are 


Fic. 3. Graph showing the region of interest in the 
T,p.4 plane for the sample problem. 
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Fic. 1. Heat generation function of q(T. p 4). 
Photograph of three-dimensional model. 
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graph of three-dimensional model. 


VOL. 
60 
j 


Intraparticle diffusion and conduction in porous catalysts—II 


Geo + + T = 8463 + 638-7 py — 14-377 
Doo Dio Pa + b,,T = 4.676 — 18-70 py 
Then 


pa (r) = 0-2500 + 0.0331 


u = — 38-07 

v 268-8 

A = 4-825 

p= 8-791 

— 0.006784 

B= 0.000 (by = 0) 


1 
«{ W (e) G, (r, €) de 


The calculations were carried out as suggested 
earlier and the partial pressure and temperature 
profiles are shown in Figs. 4 and 5. What is called 
the final solutions on these curves is a result of 
several successive substitutions in the integral 


and 
T (r) = 600-0 + 6.064 ane 


_ 
r 


1 
+ - W G, (7, €) de + 


1 
+ [we G,, (r, ©) de 


Temperature, Tir), °K 
a 
3 


Fic. 5. Temperature profiles within particle showing 
effect of linear approximation. 
~—— Final Solution 
First approximation. 


equations. The essential point here, however, is 
not the number of calculations made but rather 
the excellence of the results with the linear 
approximation, equations (8) and (9). Effective- 
ness factors are as follows : 


Final | 
solution 


Approvimate 
solution 


E.F. (on A) 


Fic. 4. Particle pressure profiles in particle for 


reactant A and product C (Exact solution). 


—— Final solution 
Pe Final solution 


E.F. (on C) 
T, (°K) 


0-364 
602 
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These calculations were made on a desk caleu- 
lator with the radius of the particle divided into 
twenty equal parts. The integrals were then 
computed using a simple trapezoidal rule. Some 
attempts were made to determine the effect of 
interval size but no conclusions were drawn. 
Attempts were also made to solve the problem on 
an analogue computer although these results 
were only partially successful. The equations 
tend to be unstable when proceeding from the 
surface inwards, Some success was achieved 
integrating the equations on an analogue computer 
starting at the particle centre and moving out- 
wards. However, two values must be guessed, 
the partial pressure and temperature at the centre 
of the particle, and this can be inconveniently 
tedious. 


SUMMARY AND CRITIQUE 


This paper is a generalization to complex 
systems of reactions of the previous paper on 
simple reactions. It was shown that the rank of 
the stoichiometric matrix gives the number of 
conservation equations which need integrating. 


An approximate method was devised for solving 


the equations and a technique was described 
which could be used to obtain temperature and 
partial pressure profiles as rigorously as desired, 
The approximation method may break down in 
some cases but it is recalled that for simple 
reactions it worked remarkably well. 

It must be stated that the methods used in 
this paper for numerical computations while 
satisfactory for two and possibly three reactions 
will probably not be suitable for larger numbers of 
reactions. Analogue computer solutions for three 
or more reactions may be less inconvenient but 
this remains to be shown. It has been suggested 
that grid methods, that is the transformation to 
finite difference approximations for the deriva- 
tions, be tried but there are considerable difli- 
culties with this because of the nonlinearity of 
the finite difference equations and lack of in- 
formation on convergence of the resulting iterative 
process. It would seem that some simpler 
approximation scheme is needed. 
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An investigation of laws governing the spouting of coarse particles* 


H. A. Becker 
(Received 5 May 1960) 


Abstract—The relationships between the variables involved in the spouting of coarse particles are 
considered. The dependent variables are divided into four categories: those dependent only on 
the gas and particle properties, those dependent also on the column geometry, those additionally 
dependent on the bed depth, and those moreover dependent on position within the bed. General 
correlations are presented for variables of the first three categories (the minimum spouting 
velocity and its maximum and the maximum spoutable bed depth), and for one of the fourth 
category (the frictional pressure gradient). Restricted correlations are given for other variables 
of the fourth category (the rate of agitation and the local gas velocity). 


Résumé —L ‘auteur considére les relations entre les variables qui interviennent pour maintenir 
en suspension des grosses particules. Il distingue 4 categories de variables dépendantes: Celles 
qui dépendent seulement des propriétés des gaz et des particules. Celles qui dépendent aussi 
de la forme de la colonne, celles qui dépendent en plus de la profondeur du lit, et celles dépendant 
en outre de la position a Vintérieur du lit. Des relations générales sont données pour les variables 
des 3 premiéres catégorices (vitesse min. de fluidisation et épaisseur min. et max. de lit fluidisé) 
et pour une variable de la 4éme categorie (gradient de pression de friction). Des relations res- 
trientes sont données pour les autres variables de la 4é¢me catégorie. (Vitesse dagitation et 
vitesse locale du gaz). 


Zusammenfassung. Dic Bezichungen zwischen den Variablen beim Wirbeln von groben 
Teilchen werden betrachtet. Die abhingigen Variablen werden in vier Gruppen eingeteilt : 
solche, die nur von den Eigenschaften des Gases und der Teilchen abhiingen, solche, die auch 
noch von den Siitulenabmessungen abhiingen, soleche, die zusiitzlich von der Héhe des Wirbelbettes 
abhiingen, und solche, die ausserdem von der Lage innerhalb des Bettes abhingen. Fiir die 
Variablen der ersten drei Gruppen (minimale Wirbelgeschwindigkeit, ihr Maximum sowie die 
maximale Wirbelbetthdhe) und fiir eine der vierten Gruppe (der Druckgradient infolge Reibung) 
werden allgemeine Beziehungen angegeben. Restriktive Korrelationen sind fiir die anderen 
Variablen der vierten Gruppe (Riihrgeschwindigkeit und Ortliche Gasgeschwindigkeit) angegeben. 


INTRODUCTION 


Tur srouren bed technique of agitating beds of 
coarse particles was discovered by Maruur and 
Gisuter [1]. They also made an exploratory 
study of the drying of wheat as a possible practical 
application of the new gas-selids contacting 
technique [2]. This study was continued and 
completed by the present author and SauLans 
[3]. In the course of the investigation it became 
clear that a knowledge of the general laws 
governing spouting was needed to guide the 
rational scale-up and design of dryers, and the 
research herein described was accordingly under- 
taken. The studies perforce transcended their 
primary aim and the generalizations to which 


they led are applicable to a variety of systems in 
broad ranges of the important variables, and have 
bearing as well on related aspects of aggregative 
fluidization. 


APPARATUS 


The type of column used in spouting consisted of an 
upright, circular cylinder with a conical base terminating 
in a gas inlet pine (see Fig. 1b). Geometrically this type 
of column is completely described by specifying the column 
diameter, D,, the gas inlet pipe diameter, Dy, and the 
angle included by the base cone, 8. Table 1 lists these 
dimensions for the five columns used in the present study. 


MATERIALS 


Table 2 summarizes the physical properties of the par- 
ticles spouted. Densities were determined pycnometrically 


* Contribution from the National Research Council of Canada, Prairie Regional Laboratory, Saskatoon, 
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with toluene. Particle diameters are expressed as the 
diameter, D,, of an equi-volumed sphere, calculated from 
the density and average weight of a particle. The cereal 
grains can be regarded as uniform particles ; the kernel 
size distribution of wheat has recently been described [4). 
The Ottawa sand was from a standard Tyler sieve fraction, 


Derinitions or Some VARIABLES AND 
Nores on Experimenta Mernops 


1. Bed Depth L 

The bed depth was measured (and enters the subsequent 
correlations) as the settled bed depth observed after 
strongly agitating a bed and then gradually shutting off 
the gas flow. (The gas used was air at room temperature.) 


2. Maximum spoutable bed depth (m.s.b.d.), L,, 


Maruur and Gisuier [1) found that for a given system 
of column, gas and particles there is a maximum bed 
depth which can be spouted. This maximum is quite 
sharply defined and is usually reproduceable within + 5 
per cent. 


3. Average fraction void, « 

The average fraction void of a bed under given condi- 
tions was calculated from the volume of the bed, the weight 
of material in it and the absolute density of the material. 


4. Fraction void at the point of incipient particulate 
fluidization, €mf- 

This characteristic, * static * fraction void was measured, 
and is uniquely defined, as the fraction void of a quiescent 
bed naturally consolidated by gradual reduction of 
agitation, as described above under * bed depth.’ 


5. Maximum fraction void, «,- 

This characteristic * dynamic * fraction void was meas- 
ured, and is uniquely defined, as the fraction void observed 
when the gas flow through a bed is increased to the point 
of maximum, coherent, dense-phase expansion of the bed, 
i.e. to the point of incipient aggregative fluidization. 


6. Minimum spouting Velocity, v,. 

The admission of an increasing flow of gas to a bed of 
particles resting in a column suitable for spouting causes 
no visible change until a central cone of particles with 
apex at the gas inlet is effectively suspended in the rising 
gas stream. The bed then expands slightly while the gas 
flow is increased a further 30 per cent or less, on which 
the entering gas jet becomes sufliciently intense to project 
a sharply bounded stream of particles upward through 
its core, forming a gas-solids jet, the ‘spout. The 
superficial gas velocity at this critical point is not repro- 
duceable to the extent that the * vield shear strength’ 
of the bed is variable, the more so when the particles are 
nonspherical or rough and the particle-to-colunin diameter 
ratio is high, promoting bridging. If, on the other hand, 
the gas flow rate through a spouting bed is gradually 


reduced until the spout just fails, the superficial velocity 
at that point is less subject to unpredictable constraints 
on the bed, and is a precisely determineable quantity 
which Maruur and Gisucer [1] call the * minimum 
spouting velocity.” 
7. Marimum of the minimum spouting velocity, v,,- 
The minimum spouting velocity increases with increasing 
bed depth and becomes maximum at the maximum 
spoutable bed depth. The experimental evidence shows (as 
will be seen) that this maximum is identical with the 
minimum superficial gas velocity causing onset of agyreg- 
ative (two-phase) fluidization. 


8. Minimum superficial velocity for onset of particulate 
fluidization, vps. 

This characteristic velocity is commonly defined as 
the superficial fluid velocity at which the extrapolated 
relation between frictional pressure drop and fluid velocity 
for a quiescent bed intersects the ordinate of equality 
between the drag and weight forces. The definition is made 
unique (at least for smooth walled columns and low values 
of the particle-to-column diameter ratio) by specifying 
that the quiescent bed is consolidated by gradual reduction 
of agitation. Then 


v at ( -A ext. 


%. Absolute local gas velocity, u 

The term “ absolute local gas velocity is here only 
meant to imply some direct, average measure of the local 
gus velocity in a bed, to be contrasted with an indirect 
measure such as the superficial velocity. The problem of 
obtaining such an * absolute measure in so complex a 
situation as that prevailing in an agitated bed of particles 
has obvious difficulties, an evidently satisfactory solution 
of which was found by using the following Pitot-type 
instrument: separate impact and static pressure probes 
were constructed from 4 in. stainless steel tubing and 
were mounted in lengths of heavy-walled brass tubing. 
The brass tubes were rigidly coupled, bringing the impact 
and static openings into apposition about 2 in. apart. 
This instrument was mounted in suitable carriers to 
probe any position in the O5ft and 0-75 ft diameter 
columns. Only beds of wheat were investigated, but 
measurements could as easily be made in beds of larger 
particles. With smaller particles, the differential pressures 
developed would be rather low, besides which the impact 
probe Reynolds number might be too small for accurate 
measurements (it is well known that in a free gas stream 
this Reynolds number must be held above 100 if the probe 
is to be used without calibration for viscous effects). The 
diameter of the impact tube for the present work was 
chosen to be comparable with the particle diameter, (but 
slightly smaller, to prevent clogging). Larger probes could 
be used, if screened, but might upset the gas and particle 
flow patterns excessively. 

Respecting successful use of this instrument, a spouted 
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bed has the saving grace that all the particles are in contin- 
uous and rapid movement. The fluctuations sensed as 
changing configurations of particles pass the probe tip 
are significantly averaged, owing to the inertia of the 
mict eter used to sense the differential pressures, 
and a constant reading is observed. The situation may 
not be much worse than when a Pitot tube is used to 
measure average velocities in a stream of high turbulence 
intensity. A wholly correct interpretation of such measure- 
ments cannot be hoped for, but interpretation in a relative 
sense can nevertheless be undertaken without qualms, 
since the nonideality of the measurements probably 
varies little. 


Resurs 
1. The four categories of dependent variables 


Maruur and Gisuier [1, 2] have described the 
general characteristics of the spouted bed 
technique of agitating beds of particles. The 
problem of quantitatively predicting the behaviour 
of spouted beds will therefore be broached at 
once, without further preamble. What follows 
is, in essence, a development and presentation of 
correlations for predicting the more important of 
the dependent variables involved in spouted 
systems. These variables are formally separable 
into four categories : 


(i) Those dependent only on the properties 
of the gas and particles: the maximum of the 
minimum spouting velocity, v,,, and the maximum 
fraction void, «,, ; 

(ii) Those also dependent on the column 
geometry: the maximum spoutable bed depth, 
L,; 

(iii) Those additionally dependent on the bed 
depth: the minimum spouting velocity, v,, the 
integral frictional pressure drop, (—A P),, and 
the maximum duration of the bed turnover 
eyele, 0, ; 

(iv) Those moreover dependent on position 
within the bed ; the frictional pressure gradient, 
(— dP dz),, the particle velocity at the wall, 
u,, the local gas velocity, u, etc. It will be seen 
that this classification delineates the most 
profitable order of consideration, that discussion 
of the variables in each category presupposes 
knowledge of those preceding it. Other approaches 
are possible but appear likely to obscure some 
basic features of spouted systems. 


2. The maximum of the minimum spouting velocity 


2.0 Introductory note. The maximum of the 
minimum spouting velocity (by Section 1) belongs 
in the first category of dependent variables, being 
a function only of the properties of the gas and 
the particles. The reason for this characteristic 
is easily appreciated on examining the relation 
between the drag forces, the weight forces, the 
reduced superficial gas velocity v/v,, and the 
reduced bed depth LL,,. typified by Fig. 1 (a). 
It is seen that at the maximum spoutable bed 
depth (L_L,, = 1) the minimum spouting velocity 
slightly exceeds the velocity (v,,y) equalizing the 
drag and weight forces in a quiescent bed, thus 


Table 1. Characteristics of columns 
Material of D,, 
construction (ft.) 
Black iron | 0-50 


Aluminium 


Glass 


* This arrangement was used in studying the onset of 
fluidization. A calming section preceded the support 
sereen holding the bed. 


indicating the existence of a terminal condition. 
Indeed, the experimental evidence (more of 
which is produced in the following Sections) 
uniformly indicates that the maximum of the 
minimum spouting velocity is identical with the 
minimum superficial gas velocity causing onset 
of aggregative fluidization. 

2.1 Direct correlation of the terminal drag 
coefficient. The aforementioned characteristics of 
the maximum of the minimum spouting velocity 
indicate that the total drag coefficient based upon 
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Fic. Ma). A typical relation be- 
tween the reduced drag force, the 
reduced gas velocity and the reduced 
bed depth for a bed of particles con- 
tained in a column shaped for spout- 
ing. This constitutes a “ phase 
diagram ~ for such a bed. Data on 
wheat in the 0-5 ft dia. glass column 
with 8 = 180°. 


equality between the drag and weight forces, 
Cy ~ + Do — p) Spu,,?. and the corresponding 
Reynolds number, Re, — D,v,, a. should be 
be uniquely related for particles of similar shape, 
(the choice of a characteristic particle dimension 
in these definitions is arbitrary, but the diameter, 
D.. of an equi-volumed sphere is convenient 
and unambiguous). The limited far 
obtained are in reasonable agreement thereon ; 
Fig. 2 indicates that the loci of data on similarily- 


data so 


shaped particles are smooth curves whose dis- 
placement with increasing nonsphericity of 
shape is in the expected direction of an increasing 
drag coetlicient. The data on spheres are closely 
represented by the equation 


Cp (1) 


The (broken) curves for the nonspherical particles 
are based on the assumption (good for moderate 


Fic. 1(b). Sketches of beds in some 

of the states represented in Fig. 1 (a): 

(1) A quiescent bed resting in a col- 
umn shaped for spouting. 

(2) A spouting bed at the point of 
minimum spouting; ie. at the 
point of extinction of the spout, 
at which v = 

A spouting bed at a higher gas 
velocity ; v > 

A slugging bed showing the 
passage of massive pockets of 
gas in a pattern unproductive 
of particle turnover ; v > v,. 


fy Siuggng 
2 / 
96 
= 
fa 
a 
Gas 
(2) (3) 
248 


An investigation of laws governing the spouting of coarse particles 


Fic. 2. A correlation of the maximum 
of the minimum spouting velocity. 

Present datae BAV © 

Maruur and O A O 

Spheres* @ O Flax seed & 

Wheat Barley © 

Corn Oats 

Sunflower seeds - 

Coffee beans 

Brucite, Shale, Gravel | 

* Spheres: peas, rape seed, mustard 

seed, Ottawa sand. 


ranges of Reynolds number) that the fractional 
excursion of the ordinate (Cp) with change 
in particle shape is independent of the Reynolds 
number, and is therefore representable in equation 
(1) by inserting an experimental shape factor 

22 + 20" (2) 

Re,, 

Appropriate values of this shape factor are given 
as the parameter in Fig. 2. 


(¢) 


Table 2. Characteristics of particles 


0-0227 
O-0218 
O-0071 
0-0053 
O-0025 


O75 


Corn 

Peas 

Barley 
Wheat 

Flax seed 
Rape seed 
Ottawa sand 


0-98 


* The densities of the cereal grains are quoted for a 
moisture content of about 0-10 Ib Ib, the natural moisture 
content at the time of the experiments. 


2.2 Correlations based on drag equations for 
quiescent beds. Since the maximum of the 
minimum spouting velocity is evidently identical 
with the minimum superficial gas velocity causing 
onset of aggregative fluidization, it may be 


anticipated that it can be predicted from the 
drag equations for quiescent beds, (the term 
“quiescent bed” is in this paper used synony- 
mously with the descriptively inferior, but 
commoner, terms “packed bed” and “ fixed 
bed"). Here the fraction void should be taken 
as the “dynamic” maximum fraction void, 
€»: and the drag forces equated with the weight 
forces, 


(— AP), = gL, (p, — p) (1 — (3) 


When this was done the above expectation was 
indeed precisely satistied by the nonspherical 
particles. The spheres (peas and rape-seed), 
however, failed in the expanded (bed) state to 
follow their drag equation for the quiescent state, 
due to a change in their fraction void function 
with the onset of bed expansion. This indivi- 
duality of behaviour is immaterial in the present 
context, since equation (1) predicts v,, for spheres 
very adequately. The present concern is therefore 
satistied by the fact that v,, can be accurately 
predicted for nonspherical particles, at least, 
given a reasonable value of «,,. adequate informa- 
tion on particle shape and a reliable drag equation 
for the quiescent state. 

The drag equations for quiescent beds can, 
however, be used for the same end in a slightly 
different way, such that the spherical particles 
also become obedient to a general rule. The 
minimum superficial gas velocity for the onset 
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of particulate fluidization, v,,, can always be 
calculated from the drag equations given reason- 
able values of the appropriate fraction void, 
«ay: It is therefore of interest to inquire whether 
v»y and v,, are not directly related, perhaps by a 
simple proportionality. Table 3 indicates that 
this ts indeed the case and suggests the simple 
rule that 

Um 1:25 (4) 


The Table also gives the experimental values of 


the characteristic fractions void, and 


Table 3. Critical data on the agitation 
of particles 


Um 


1260 O95 
123 
1-33 — | - 
126 | 0-380 
122 0-371 
110 | 


112 O-395 
| 
* The effective fraction void of a barley bed cannot be 
measured, owing to the peculiar structure of the barley 
kernel. 


A rather exhaustive investigation of drag in 
quiescent beds of spheres, wheat and corn was 
made in the course of the present study, and it 
may be of interest to record the results here. 
The data on spheres (« — 0.35-0-44) are very well 
correlated by the equation (for D, D, < 0-05) 
K &exp[—42(1 — 24Re 4 186? Re® (5) 
The data on wheat (¢€= 0-37-0438) and corn 
(« = 0-37-0-46) are well correlated by the 
equation 
K @exp[— 42(1 Cy Re 4+ C,Re®, (6) 
as may be seen in Fig. 3, but can be more simply 
represented by the restricted equations 
Keexp[ #2(1 64 

60 Re < 250 
for wheat, and 


K &exp [— 42 (1 — = 38 
150 Re < 1000 (8) 
for corn. 

The foregoing equations will not be justified 
here ; the investigation of drag in quiescent beds 
is being continued and will be the subject of 
another article. 


Fre. 3. A correlation of the maximum of the minimum spouting velocity. 


| 
Cp K/Re® xX 


Spheres 4:5. 5 Re | Cpe ®exp[—42(1 — 


Non-spheres | Re 


Cp — (420 — 


VOL, 
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3. The marimum spoutable bed depth 


3.0 Introductory note. The maximum spout- 
able bed depth (m.s.b.d.) comprises the second 
category of dependent variables, according to 
Section 1; in contrast with the maximum of 
the minimum spouting velocity, the m.s.b.d. is 
dependent on the dimensions of the column 
used in spouting as well as on the properties of 
the gas and particles. Quantification of the 
attributes of the m.s.b.d. thus introduces the 
novel and interesting problem of the modelling 
of spouted beds. 


3.1 Phenomenological note. The m.s.b.d. marks 
the closure of the stability of spouting with 
respect to changes in the gas flow rate. At 
lesser bed depths increases in the gas velocity 
above the minimum required for spouting aug- 
nent the rate of bed turnover and cause the 
spout to mount higher above the bed, (see Fig. 1b). 
Ultimately, at an instability point characteristic 
of the bed depth they cause supersession of 
spouting by less organized agitational patterns 
such as slugging (with coarse particles) and 
aggregative fluidization (with fine particles). At 
small bed depths the gas velocity can be increased 
almost indefinitely without causing instability. 
But, as the m.s.b.d. is approached, the loci of 
the points of instability and minimum spouting 
draw together (Fig. la) and finally meet at the 
m.s.b.d. itself. The spouting bed may then be 
said to be at a critical phase transition point, as 
illustrated by Marner and Gisuters’ * phase 
diagrams’ [1] and by Fig. 1 (a). 

The m.s.b.d. (enlarging on remarks in the 
preceding Section) also marks the point of 
maximum fractional expansion of beds, both prior 
to and during spouting. The physical appearance 
of a spouting bed is, to a first approximation, 
independent of the bed depth, and dependent 
only on the reduced level, z L,,, at which the 
bed is observed. At low reduced levels the bed is 
in a compact state, the upward flow of gas being 
largely accommodated by the spout which is 
describable here as a high-velocity jet piercing the 
bed and supporting a thin population of particles. 
At values of z/L,, of 0-3-0-4 there is a noticeable 
loosening or expansion of the annulus of the bed. 


This loosening grows with increasing value of 
2/L,,, and becomes maximum at z= The 
condition of the bed at levels near z = L,, is, 
in fact, that of incipient aggregative fluidization, 
the uppermost levels of the bed having reached 
the limit of coherent, dense-phase expansion, 
It is therefore evident that the gas flow must 
be nearly uniform at z — L,, and that the super- 
ficial gas velocity must be practically identical 
with the minimum velocity for onset of aggrega- 
tive thuidization. 

The m.s.b.d. therefore characterized by 
maximum gas loading of the annulus, by maximum 
solids loading of the spout, and by consequent 
degeneration of the agitational pattern into aggre- 
gative fluidization or slugging. 

The foregoing deductions are further borne 
out by gas velocity profiles presented in Section 7. 
That the fractions void of spout and annulus 
approach each other at z — L,, was clearly visible 
in high-speed photographs taken of two-dimen- 
sional spouting in a narrow, rectangular column. 

3.2. A hypothesis for correlating the maximum 
spoutable bed depth. The m.s.b.d. can, in principle, 
be related to the independent variables through 
force and energy balances. However, an exact 
statement of any kind would involve knowledge 
of momentum and energy transformations in the 
spout difficult or impossible to acquire experi- 
mentally. One consequently seeks a simplification 
which retains the useful skelton of a force or 
energy balance but which excludes those elements 
containing dependent variables of higher order 
than the m.s.b.d. itself, (i.e., variables of the third 
and fourth categories in Section 1). The simpli- 
fication used in this paper is formulated as follows. 

It is noted that the formation of a spout is 
accompanied by a significant decrease in the 
frictional pressure drop through a bed. (At the 
m.s.b.d. the frictional pressure drop through a 
quiescent bed on the verge of spouting essentially 
equals the weight stress at the foot of the bed, 
Equation (3), but when the spout breaks through, 
the decrease in gas flow through the annulus 
coincident to partial venting of the total gas flow 
through the spout causes a notable fall in the 
pressure drop, Fig. la). Thus the notion of a 
force balance suggests the hypothesis that 
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change in gas 
momentum in (9) 
the bed 


decrease in integral 
pressure drop at the> o 
advent of spouting 


This statement is at least qualitatively correct, 
and is of a type facilitating analysis of a situation 
sulliciently complex that unguided dimensional 
analysis is so ambiguous as to be nugatory. 

The hypothesis may be written symbolically 
in the form (with reference solely to beds of 
maximum spoutable depth) 


PQ, (Um — 4) (10) 


The velocity of exit, u,, is the momentum-mean 
velocity within the bed, just prior to exit, and will 


for this reason be expressed as u, = &, v,,. where 


k, is the factor times which the momentum-mean 
velocity exceeds the superficial velocity within 
the bed. 
putting Q, = Deu, 4 and = »v,, 
one obtains 


Thus, rearranging equation (10) and 
D2 DZ 


1 —k,(D, DP 


(l —«,,).@ (11) 


This rather unwieldy expression can be relaxed 
by taking account of the following facts: (i) the 
present experimental data indicate that the 
reduced pressure decrement, ®, is proportional 
to powers of the ratios Dy and ; (ii) the 
drag coeflicient C, is given in terms of particle 
shape and Reynolds number by equation (2) ; 
and (iii) the fraction unvoid (1 — «,,) is uni- 
quely determined by particle shape and (possibly) 
Reynolds number, Hence equation (11) can be 
written, allowing for residual effects of the flow 
regime and particle shape, 


3.3 Correlation of 
maximum spoutable bed depths. Examination of 
the abundant experimental data on wheat, 
(Fig. 4). indicates that the exponent m in equation 
(12) has the value 1-76. Let it be assumed (as 
justified by the results) that this value of m is 


4 


[f 


experimental 


(Re,,) (12) 


data on 


universal, Figs. 5(a) and 5(b) then indicate 
that the exponent n is expressible as 


n= Léexp(— 0.0072 Re,,) 


Present dato 


Mathur and Gishier 


Fic. 4. A correlation of the maximum spoutable bed 
depth for wheat beds. 

Ging Vin. B Pin. O@ 

The final correlation, Fig. 5(c), then follows 

directly on determining the optimum exponent 

of the shape function, f(¢). The equation of the 

correlating curve is 


Re, 42 (14) 


The correlation is satisfactory for the range of 
particle shapes and Reynolds numbers so far 
investigated. Higher Reynolds numbers are 
probably uninteresting in practice, since extremely 
coarse particles would be involved. The limiting 
behaviour at low Reynolds numbers cannot be 
forecast with any accuracy, and data extending to 
Reynolds numbers as low as 0-5 are needed to 
fill this gap. 

3.4 A correlation for low Reynolds numbers. 
An alternative interpretation of the present data, 
(Fig. 6), indicates a simplified behaviour at 
Reynolds numbers Re,, smaller than about 100 
and leads to the correlating equation for spheres 


44 16 10 
Re, 


This equation is, to be sure, incompatible with 


(15) 
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Re,, 
100 


N 


Res 670 


Sona 


Re, 


Fic. 3(a). Correlation of the maximum spoutable bed depth for corn, 
rape-seed and Ottawa sand. Data on the 0-5 ft dia. columns. 
Fic. 5(b). Correlation of the exponent n. 

Fic. 3(c). Generalized correlation of the maximum spoutable bed 
depth. For key to data see Fig. 2. 


equation (14), and fails emphatically at Reynolds 
numbers above 100. Nevertheless, it gives an 
undeniably fair portrayal of the low Reynolds 


Le 
(=, 


Present doto 
Mathur and Gishier 


Fic. 6. A correlation of the maximum spoutable bed 
depth for Reynolds numbers Re,, smaller than about 
100. 


number behaviour and strongly suggests that the 
m.s.b.d. varies in direct proportion to the gas 
inlet particle Reynolds number Re, when the 
inlet flow regime is viscous with respect to the 
particles. 

3.5 The significance of the geometrical group Q. 
It can be plausibly argued that the group 
= (L,, D,)(D, which appears in all 
the correlations, is a geometrical factor charac- 
terizing the ratio of the axial and radial resistances 
to gas flow in the spread of the gas from the spout 
through the annulus. This explanation has the 
necessary merit that it does not overtly contradict 
the constancy of the exponent, 1-76, on the ratio 
D,. 

3.6 Limiting effects of the inlet-to-column diameter 
ratio. The group 1 — k, (D, D,)? in equation (11), 
eliminated from the foregoing analysis, has still 
an important significance with respect to the 
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stability of spouting. Clearly, if k, is larger than 
unity this group will become zero at a value of 
D, D, smaller than unity. In facet, evidence 
presented in Section 7.1 indicates that the value 
of k, in a turbulent regime is approximately 
eight, giving a critical value of D, D. of 0.35, 
This means that beds will not spout if Dy D, 
exceeds 0.35; at — 0.35 the m.s.b.d. is 
zero. A fairly accurate physical interpretation 
of this prediction is that with D, D, greater than 
0.35 the gas velocity in the bed exceeds the 
minimum velocity for onset of aggregative 
ausnuzation before the inlet velocity exceeds the 
free-fall velocity of the particles. 

3.7 Note on a Mathur and Gishler correlation. 
It is of interest to note that Marnun and Gisuien 
[1] have proposed a correlation for the minimum 
spouting velocity which should yield an expression 
for the m.s.b.d. if one substitues L — L,, and 
Uy Um: The ircorrelating equation gives, after 
rearrangement (see p. 164 of their article), 


2 2a 
| Fal Fal Cp = constant (16) 


This equation contains the basic elements of 
equation (14). The “ constant on the right is, 


7. 


however, a function of column diameter (private 
communication). 


4. The minimum spouting velocity 


4.0 Introduction. The third category of depen- 
dent variables, by Section 1, comprises those, 
such as the minimum spouting velocity and the 
integral frictional pressure drop, which are func- 
tions of the bed depth as well as of the column 
geometry and the properties of the gas and the 
particles. However, only the minimum spouting 
velocity will be discussed at this point as the 
integral pressure drop most conveniently 
interpreted with the frictional pressure gradient 
in the fourth category of variables. The following 
discussion introduces a concept which ts basic 
to the modelling of spouted beds. 

4.1) The similarity principle of spouted beds. 
Investigation of variables of the third and fourth 
categories produced clear perception of the 
following similarity principle : 

Beds of given particles spouted with a given 
gas in geometrically similar columns behave 
similarly when the reduced bed depths LL, are 
the same. 

4.2. The equation of the minimum spouting 


A correlation of the minimum spouting velocity 


for wheat spouted in geometrically similar columns. 


Dy/D, = 0-172 
D, = 240 ft 


D, = 1-00 ft 
D, = 0-75 ft 
D, = 0-50 ft 


B= 108° @ 
B 
B= 
B=120° 
B= 120° 
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Fig. 8. 


A correlation of the minimum spouting velocity for 


wheat spouted in the 0-5 ft dia. glass column with 8 = 180°. 


velocity. In the particular case of the minimum 
spouting velocity the above similarity principle 
predicts that for a given particle and gas and 
geometrically similar columns 


Um D,,) 
The validity of this relation is illustrated by Fig. 7, 
which shows a correlation of data on wheat 


(the m.s.b.d. in the 2 ft diameter column was too 
large to be determined experimentally and was 


(17) 


therefore calculated from equation 14). The 
relation appears to be unaffected by variations 
in the base-cone angle between 90° and 180° 
(the m.s.b.d. was also unaffected by these 
variations). 

Fig. 7 suggests the explicit form of equation (17) 


L 
1 4) 
8s n(; 


“m 


2 (18) 


Um 


This relation was universally exhibited by the 
present data, excepting those for very low bed 


eol- 
‘ 0065,0085 


20- 
20-15 


40-085 


= ~ 
0-085e 


y= 150 


| 


For key to data see Fig. 2. 
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Re” 298. [+ 5 


Fic. 10. 


depths ; the lower limit of validity seems to be 
about L D, — 1. Thus the problem of correlating 
the minimum spouting velocity reduces to 
determining the behaviour of the numerical 
coeflicient s in equation (18), 

4.3 Correlation of the coefficient, s. The coetli- 
cient s decreased with increasing value of the 
geometrical ratio D, D,, as may be seen in Fig. 8. 
Examination of the present data showed that 
the product sD, D, is exactly constant for a 
given particle excepting certain of the data for 
low Reynolds numbers. Variations in this 
product depend on the Reynolds number Re,, 
(Fig. 9a) at least at Reynolds numbers above 100. 
The correlation is improved somewhat by postula- 
ting a 2 3 power dependence on the shape factor 
S(¢), Fig. 10b. Fig. 10a indicates that the gas 
inlet particle Reynolds number Rey is a possible 
criterion for deviation at low Reynolds numbers, 
but Fig. 10(b) suggests that the inlet pipe 
Reynolds number may be a better one; the 
data for Dy typ > 30000 are uniformly well 
correlated by the equation 


s 0.0071 (7) [f Re (19) 


A general correlation of the coeflicient s. 


For key to data see Fig. 2. 


The validity of the Reynolds number criterion 
Dy tig p should, however, be further investigated 
at Reynolds numbers Re,, smaller than those in 
the present study. 

Combination of equations (18) and (19) gives 
as the equation of the minimum spouting velocity 


1 + 0.0071 


[ f(¢) Re,,°2 In ) (20) 
Figs. 9 and 10(b) indicate that this equation is 
quite generally valid at Re,, > 100, while at 
10 < Re,, << 100 it is valid only if Dy uy p 
30000. However, since no firm justification can, 
on the basis of the data so far obtained, be 
offered for the criterion Dy tgp perhaps 
Fig. 9(b) should be interpreted in the following, 
simplest sense : Equation (19) is valid in the range 
10 < Re, < 100 if Dy D, < 0-1. 

4.4 A correlation for very low bed depths. 
The minimum spouting velocity seems to become 
independent of the properties of the particles as 
LD, falls below about unity. The following 
equation correlates data on all the particles as 


i 
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spouted in the 05ft diameter glass column 
with 8 — 180 : 


1-2 ( exp | (21) 
Um dD. dD, 
Other data indicated that increasing the cone 
angle to 90° increases the numerical coeflicient 
of the r.h.s. slightly. 

There should be little practical interest in 
beds with LD. <— 1, for at such low bed depths 
agitation is feeble and of poor quality and the 
greater part of the gas flow is localized in the 


region of the spout. 


5. The frictional pressure gradient 


5.0 Note on variables of the fourth category. The 
remaining Sections consider variables of the 
fourth category, i.c. variables which are not only 
functions of the gas and particle properties, the 
column geometry and the bed depth, but also of 
position within the bed. In the cases of the 
frictional pressure gradient and the particle 
velocity at the wall position in the bed is ade- 
quately specified by stating the level in the bed, 
expressed by the reduced height 2 L,, above the 
gas inlet. But for the local gas velocity radial 
variations are not negligible and one must also 
specify the reduced radial position 7 r,. 

5.1 A hypothesis for the gas flow in the annulus. 
The frictional pressure gradient reflects the local 
gas flow in the annulus, and is determined by the 
drag equations for quiescent beds when the 
fraction void and superficial gas velocity in the 
annulus are known as a function of bed level. 
The present data show that the frictional pressure 
gradient is satisfactorily accounted for by the 
following simple hypothesis as to the gas flow 
in the annulus : 

The mean superficial gas velocity at any level 
in the annulus is approximated by the equation 
of the minimum spouting velocity (equation 20) on 
substituting v, for v, and = for L: 

Va 
1 + 0.0071 x 


[f (6) In (22) 


5.2 A simplified equation for the frictional pres- 


sure gradient. With the above nypothesis the 
equation can be written, linking the gas velocity 
in the annulus with the frictional pressure 


gradient, 
b 


1 (28) 


g(p, —p)(l —«) dz 
where « is an average value of the fraction void 
(say, [<p +e ] 2) and the exponent b is given 
by the drag equations for quiescent beds ; e.g. 
by equation (6), b — 1-65 for wheat. Equation 
(23) is approximate, not exact, since the effects 
of gradients in the fraction void and in the super- 
ficial gas velocity are neglected. 

Nevertheless, it provides reasonable estimates 
of the frictional pressure gradient when equation 
(22) is substituted for v, v,,. 

5.3 The integral frictional pressure drop. 
The integral frictional pressure drop (a variable 
of the third category) is easily calculated from 
the frictional pressure gradient and, as noted 
earlier, is always somewhat smaller than the 
weight of the bed on unit area. A_ practically 
significant fact, however, is that a pressure drop 
equal to weight of the bed on unit area must be 
allowed for to start a bed spouting, this being the 
maximum pressure drop attained prior to forma- 
tion of the spout. 


6. The rate and quality of agitation 


6.0 Introductory note. An idealized definition 
of spouting might specify a short, uniform bed 
turnover cycle consisting of downward flow of 
the annulus, capture of particles by the spout 
at the base of the bed, nearly instantaneous 
upward projection of these particles in the 
spout and their fall back into the annulus at the 
top of the bed. In actuality, however, particles 
enter the spout at all levels in a bed, and the rate 
of turnover therefore decreases from top to bottom 
of the bed. This may be an adverse factor under 
some conditions, and it is well to know the 
conditions which favour approach to ideality. 

Bed turnover can be studied in two simple 
ways: (i) the particle velocity at the wall may be 
directly observed, and (ii) a marked pellet 
(or pellets) may be dropped into the bed adjacent 
to the wall and the lapse of time to its 
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reappearance at the top of the bed noted. The 
particle velocity at the wall measures the velocity 
of the annulus there (radial variations in the 
annular particle velocity should be small, since 
they require shearing of the annulus), and the 
decrease in annular velocity with decreasing 
height above the gas inlet measures the local 
flow of particles into the spout. The tracer 
method, on the other hand, measures the 
maximum evele time in the bed (a variable of 
the third category). Both these methods were 
used in the present work to study turnover in 
beds of wheat. Further studies are needed to 
determine the generality (and possible moditica- 
tions) of the following correlations. 

6.1 Correlation of the particle velocity at the wall, 
The particle velocity at the wall of wheat beds 
spouted in the 0.5 ft diameter glass column was 
correlated using the similarity principle of 
Section 4.1 and the hypothesis as te gas flow 
in the annulus in conjunction with plausible 


Fig. 11. A correlation of the particle velocity at the 
well for wheat spouted in the 0-5 ft dia. glass column 
with 8 = 180° and Dy D, = 0065-0285. 


dimensionless groupings of the variables (Fig. 11). 
The correlating equation is 


where B= 0.055 and at 2 L, > 0-25, 
and B 0-21 and b = 2 at L,, > 0-25. Sub- 
stitution of equation (22) for uv, v,, gives 
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5/2 


+ 0.0071 (7) (25) 


for: L,, > 0-25. Fors L,, < 0-25, the numerical 
coetlicient of the r.h.s. is 0-142, and the exponent 
of 2 L,, is 5 2. 

Equation (25) should provide a safe basis for 
gauging the effects of the independent variables 
on the agitation of any given particles in a given 
column. It may not, however, predict fully 
the effects of the Reynolds number Re,, and the 
particle-to-column diameter ratio D, D,. 

6.2. Correlation of the maximum turnover cycle 
time. The maximum period of the bed turnover 
eyele, @,, was studied in wheat beds 2-6 ft deep 
spouted in all the columns listed in Table 1, 
with D, D, 0.17 and 90 and 180°. The 
data are satisfactorily correlated by the equation 


0, Wy (hg) (Sn) (26) 


This equation was developed merely to provide 
a basis for treating certain problems encountered 
in the design of spouted wheat dryers, but 
probably contains the germ of a more general 
correlation. The maximum cycle time can be 
computed from the particle velocity at the wall 
but difficulties are encountered at the base of 
the bed, where the movement is chiefly inward 
and the downward velocity approaches zero 
However, it can be shown that equations (25) 
and (26) are roughly compatible. 

6.3 The effect of gas flow rate on agitation. The 
foregoing correlations are for the condition of 
minimum spouting. Increasing the gas flow rate 
above the minimum increases the rate of agitation 
roughly in proportion to the gas flow ratio vu v,. 


7. Local gas velocity patterns 

7.0. Introductory note. Although the basic 
gas flow pattern in a spouted bed is intuitively 
obvious, it is still of interest to study its details 
quantitatively. Basic data may be obtained for 
analysing the spread of the gas into the annulus, the 
local effectiveness of gas-solids contact may be 
judged and subtle effects of changes in the 
operating variables may be observed. 
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7.1. Some observations on the measured local 
velocities. The experimental data on wheat 
(no other particles were studied) indicate that the 
measured local velocity in a spouted bed at 
levels approaching = — L,, is of the same magni- 
tude as the terminal free-fall velocity of the 
particles in an infinite fluid. Hence local velocities 
are appropriately expressed in terms of the ratio 
uu, Remembering that bed levels near = — L,, 
are in a state of incipient aggregative {luidization 
this fact is striking. The identical situation must 
exist in the dense phase of thuidized beds. How- 
ever, the finding may well prove to be valid only 
for non-viscous flows. 

For the present case (of non-viscous tlow) it 
is of interest to compare the terminal, free-fall 
velocity with the maximum of the minimum 
spouting velocity and to determine the compat- 
ibility of their measured ratio with several 
hypotheses as to the effective local value of 
the gas velocity in a quiescent bed. The terminal 
velocity of the wheat was found to be 24 ft sec 
(caleulated from the author's correlation [5] after 
experimentally checking the shape factor by 
free-fall experiments in aleohol and water) and 
the maximum of the minimum spouting velocity 
was 3-1 ft sec, giving a ratio u, v,, of eight. If 
volume effects alone are considered, the effective 
local velocity should be 1 «,~5 2 times the 
superticial velocity. If, however, the local velocity 
is determined by the minimum free-flow area 
between contiguous spheres, the maximum velo- 
city ratio is y3[\ (3) — (7 2)]— 11. The meas- 
ured local velocities are thus seen to be of the 
correct order of magnitude, the second hypothesis 
being in best accord with experiment. 

In view of the foregoing, it can be plausibly 
argued that the measured local velocities closely 
approximate both the momentum-mean and 
residence time mean velocities in a turbulent 
regime. Certainly much of the void volume of 
a bed must be nearly stagnant with respect to 
gas tlow in a turbulent regime, ic. in the sense 
that fluid is caught in wake regions of low velocity. 
The assumption sometimes made that the local 
velocity approximates the void space-mean 
velocity appears to be erroneous. 

7.2 The axial velocity decay. Fig. 12 shows 


data on the decay of the axial velocity with 
increasing height above the gas inlet in beds of 
near-maximum spoutable depth. At these bed 
depths the data fall on a curve characteristic of 
the m.s.b.d.; at lower bed depths the initial 


Fig. 12. Decay of the axial gas velocity in spouted 
wheat beds. 
vine L/Ly =10 
6in OF O8e 


76e 


Fig. 13. Decay of the axial gas velocity excess in 
spouted wheat beds. 
D, =%in @ L/Ly = 10 2 
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portion of the curve is depressed in proportion 
to the decrease in inlet velocity. The correlating 
equation for bed levels in the range 0-1 < 2 L,. < 
0-45 is 


(27) 


and is (but not because of any very strong 
Huid-dynamical analogy) similar to the decay 
law [6] of a free jet 


u 
(28) 


D, 


No 


Fig. 13 shows the form of the approach to an 
uniform velocity distribution at exit ; the decay 
law of the axial velocity excess is exponential : 


u 
(29) 


1 746 exp ( 


7.8 Radial velocity profiles. Fig. 14 shows 
radial velocity profiles for beds of near-maximum 
spoutable depth spouted in the 0.5 ft diameter 
column, The behaviour exhibited is unremarkable 
excepting the development of a pronounced 
erest at the wall. Data on the 0.75 ft. diameter 
column (Fig. 15) however, show an unexpected 
complexity : each profile has an undulating shape 
ending at the higher bed levels with the same 
pronounced crest at the wall. When a series of 
closely-spaced protiles of this type are placed 
one above the other, the evolution of the standing 


wave forms is easily traced through the bed. 
There is no readily apparent explanation for this 
phenomenon which was completely absent in the 
0.5 ft diameter column. Possibly it is something 
that is favoured by small values of the particle- 
to-column diameter ratio. 


0-8 


Fic. 14. Radial gas velocity profiles for wheat beds 
spouted in the O-5 ft dia. column with 8 = 180, 


Velocity profiles for intermediate bed depths 
are nearly superimposable on those in Figs. 14 
and 15 but the trough between the wall and the 
spout tends to be more pronounced. At low bed 
depths the gas velocity at the wall is markedly 
depressed, 


oe 0 


Fig. 15. Radial gas velocity profiles for wheat beds spouted in the 0-75 ft dia. column with B = 120°. 
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7.4 Effects of the gas flow rate. Data presented 
in the foregoing are for beds at the point of 
minimum spouting. The effect of increasing the 
gas velocity above the minimum is to increase 
the gas flow through the spout without signific- 
antly affecting the flow through the annulus 
(the constancy of the frictional pressure drop 
points to the same conclusion). The excess. gas 
flow is therefore ineffectively contacted with the 
particles, and serves only to augment the rate 
of agitation. 


CONCLUSION 


The foregoing account is, hopefully, sufliciently 
complete to expedite design calculations for 
almost any feasible application of spouting that 
may arise in practice. Gas and particle properties 
will usually be given. The scale of the application, 
consideration of the qualities of agitation and 
gas-solids contact and an economic balance will 
usually determine the size of the column and 
the bed depth (the use of batteries of columns 
may be considered in place of one large one). 
The bed depth determines the maximum gas 
delivery pressure and the minimum spouting 
velocity. In practice, a gas flow up to 30 per cent 
in excess of the minimum should be available 
to ease the starting of agitation. The use of a 
90° base cone will normally obviate stagnation 
of the bed around the gas inlet. 

Spouting is, as suggested by Maruur and 
GisuLer [1], primarily a technique for agitating 
particles too coarse to fluidize well. Indeed, there 
may be not competition with fluidization at all 
since it has not yet been shown that particles fine 
enough to fluidize really well will truly spout. 
It has been seen that spouting takes place at 
superficial gas velocities of the order of the 
minimum for the onset of aggregative fluidization 
and becomes unstable at velocities of the same 
order. Thus, with fine particles spouting would 
anyway yield to aggregative fluidization at all 
but perhaps uninterestingly low gas velocities. 
However, the use of a constricted gas inlet may 
still have a beneficial effect on fluidization per se. 
This does not preclude the manifestation of 
phenomena resembling spouting in beds of fine 
particles ; these should simply not be confused 


with genuine spouting whose characteristics have 
been quantified in the preceding pages. (Note : 
Strongly oblate or prolate bodies have a marked 
tendency to channel. If such a tendency is 
present, a pseudo-spout (an internal channel) 
may form under any conditions. Such pseudo- 
spouts, unlike real spouts, have a_ peculiar, 
monotonous stability, e.g. they are insensitive 
to the inlet-to-column diameter ratio and are 
markedly insensitive to bed depth, showing no 
sharply defined “ m.s.b.d.”, tending merely to 
drift to the wall at some large, variable bed 
depth. This behaviour was actually observed in 
some beds of barley and flaxseed. Pseudo- 
spouting is a rather feeble form of agitation ; the 
pseudo-spout is frequently nothing more than 
an empty vent owing its stability to the bridging 
tendency of the particles). 

The concept of maximum spoutable bed depth 
provokes speculation as to the existence of a 
maximum fluidizable bed depth for aggregative 
fluidization. The m.s.b.d. decreases with decreas- 
ing Reynolds number as the viscous regime is 
entered. Should not the maximum fluidizable 
bed depth decrease with increasing Reynolds 
number as the viscous regime is departed 
from ? 

Finally, to emphasize that spouting and aggre- 
gative fluidization are fundamentally different 
phenomena, consider the basic mechanisms of 
agitation: in spouting agitation is caused by an 
intense, continuous jet of gas piercing a quiescent 
bed and yielding its energy to generate a well- 
defined stream of particles, the spout, whose 
upward motion is balanced by the downward 
motion of the surrounding annulus. In aggregative 
fluidization on the other hand, agitation is largely 
due to the stirring action of trains of gas cavities 
or bubbles rising through the continuum of 
loosely contiguous particles. In the former, jet 
momentum governs; in the latter, a buoyancy 
acting on gas volumes exceeding the saturation 


gas flow. 
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Novation 
Dimensions in absolute units 


length 
time 
column (inside) diameter L 
gas inlet pipe (inside) diameter L 
particle diameter, expressed as the diameter 
of an equi-volumed sphere L 
drag force on a particle ML @ 
7 DS(— A Py — ©) in beds of part icles 
7 DAg (po, — p) 6 at equality with the weight 
forces 
local acceleration due to gravity Le 
= bed depth L 
maximum spoutable bed depth L 
pressure M Le 
frictional pressure drop M Le 
vm 4 minimum gas flow rate caus- 
ing spouting at the maximum spoutable bed 
depth Le 
radial co-ordinate originating at the column 
anis L 
D2 radius of column L 
Dy/2 © radius of gas inlet tube L 
local gas velocity in a bed Le 
Um (D, inlet gas velocity at 
Le 
momentum-mean gas velocity  L,, L @ 
terminal free-fall velocity of a particle L @ 
axial co-ordinate originating at the gas inlet L 
maximum duration of the bed turnover cycle 
0 
gas viscosity M LO 
gas density ML 
particle density ML 
superficial gas velocity Loe 
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superticial gas velocity in the annulus Lp@ 
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Dimensionless variables 


numerical exponents 


= numerical coeflicient 
= cocflicients in quadratic drag law 
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mum spoutable depth 


= experimental particle shape factor 


D2 particle sphericity (S — particle 
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Hold-up in packed liquid-liquid extraction colymns 
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Abstract—A generalized correlation for hold-up of dispersed phase in packed liquid extraction 
columns has been derived and used to correlate the data of the present investigation in 2 in. 
column along with the reported data in 3 in. to 12 in. columns for a wide variety of liquid systems 
using different sizes of Raschig rings, Berl Saddles and Lessing rings. 


Résumé—Une relation générale exprimant la rétention de la phase dispersée dans une colonne 
& garnissage extraction liquide liquide a été établie ; elle est utilisée pour comparer les données 
de la présente étude, obtenues avec une colonne de deux pouces, avec celles relatives & des colonnes 
de 3 a 12 pouces, pour une grande variété de systémes, et pour différents garnissages anneaux 
Raschig, selles Berl, et anneaux Lessing. 


Zusammenfassung bine verallgemeinerte Beziehung fiir den hold-up der verteilten Phase 
in Faillkérperkolonnen der Flissig-Flissig-Extraktion wurde abgeleitet und dazu benutzt, die 
vorliegenden Ergebnisse in ciner 2 in. Kolonne mit den Ergebnissen in 3-12 in. Kolonnen fiir 
einen weiten Bereich von fliissigen Systemen unter Verwendung von Raschigringen, Berl-Satteln 


und Lessingringen verschiedener Grosse korrelieren. 


INTRODUCTION 


THe DEVELOPMENT of a satisfactory method for 
prediction of the performance of packed extrac- 
tion columns will almost certainly involve the 
estimation of the interfacial area of the contact 
of the two phases. This may be estimated from 
the determination of the hold-up of the dispersed 
phase in the packing and of the mean droplet 
size. A study of droplet behaviour in packed 
columns has been made by Lewis ef al. [1]. 
Hold-up in packed columns unlike that in spray 
columns cannot be determined accurately from 
the measurements of pressure drop since allowance 
must be made for the net thrust upwards or down- 
wards on the packing. Very little information 
has appeared in the literature on the determina- 
tion of hold-up in packed columns, although it 
has been observed that in the eperation of liquid— 
liquid extraction columns a definite amount of 
dispersed phase is always present in the interstices 
of the packing and the amount of dispersed phase 
thus held up is dependent upon the flow of the 
two phases in the column. The earliest attempt 
on the measurement of hold-up was perhaps 


that of Arre: and Excry [2], carried out in the 
course of their work on the rate of extraction of 
benzoic acid from aqueous solutions. GAYLER 
and Prarr [3, 4] have reported some data on 
hold-up and pressure drop in packed columns. 
Data on hold-up at flooding have been correlated 
by Venkataraman [5]. More recently GayLer, 
et al. [6] measured hold-up below loading and 
correlated the characteristic velocity in terms 
of the physical properties of the liquid systems 
and the packing characteristics. The present 
investigation was aimed at developing a satis- 
factory and generalised correlation for held-up 
below flooding in packed liquid-liquid extraction 
columns. 
EXPERIMENTAL 

Apparatus 

The extraction column used in the present 
investigation was a 2in. id. pyrex glass tube, 
3 ft in length fitted with expanded end pieces, 
6 in. id. at the bottom and 4} in. id. at the top, 
described by Venkataraman [5]. A schematic 
diagram of the experimental set-up is shown in 
Fig. 1 for ready reference. 
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Packings 

The physical characteristics of unglazed Raschig 
rings. Berl saddles and Lessing rings used as 
packings are given in Table 1 


Table 1. Packing characteristics 


Outside 
Packing aie 


| cine | 


0-363 1582 
206-0 O58 
O25 O874 185-0 0-635 
O187 283-0 0-69 
& | — | 2280 !0-65 


365 


oO 194 
123 


Zin. Raschig rings 
} in. Raschig rings 
2? in. Lessing rings 
jin. Lessing rings 
din. Berl Saddles 


Liquid systems 

Hold-up determinations were made in conjunc- 
tion with the study of extraction of acetone from 
aqueous solutions by kerosene, butyl acetate and 
benzene. The aqueous phase (continuous phase) 
was always saturated with the solvent phase 
(dispersed phase). Similarly the solvents used 
were saturated with water. The average pro- 
perties of the continuous phase and discontinuous 
phase are given in Table 2 for the various liquid 
systems studied in the present investigation. 


S. Lappna 


Seporotor 


Table 2. Average physical properties of 


liquid phases at 30°C 


(dynjem) 


(g/em*) (g/em) cP cP 


Kerosene | 


acetone water 
Butyl acetate 
acetone water 
Benzene | 
acetone water 


OOS OSTE OSE 


0-656 


EXPERIMENTAL 

In order to ensure that each feed solvent was 
saturated with the other, a layer of water was 
maintained on the bottom of the solvent layer 
in the solvent feed tank and a layer of solvent on 
the top of the fluid in the aqueous feed tank. 
Feed flow of both phases was motivated by means 
of constant air pressure in the feed tanks. The 
column was filled with the continuous phase 
and the rate adjusted to the desired value as 
observed by the orifice meter in the line. The 
dispersed phase flow was next turned on and 
brought up to the desired value, the interface 
position having been stabilized in the middle of 
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Hold-up in packed-liquid liquid extraction columns 


the separating section at the top. When con- 
ditions became steady as evidenced by the 
constancy of the pressure drop readings, the 
aqueous and solvent tlow rates were stopped 
simultaneously. It was observed that the inter- 
face at the top descended slowly due to the 
accumulation of the dispersed phase held up in 
the column pushing down the continuous phase. 
When no further movement of the interface was 
observed, the continuous phase supply was 
turned on gradually and the dispersed phase 
liquid collected at the top was allowed te be 
displaced out of the column in a receiver. The 
continuous phase supply was cut off exactly 
when the interface at the top regained its original 
level. The volume of the dispersed phese collected 
was determined and the ratio of the volume of 
the dispersed phase present in the interstices of 
the packing between the interface and the packing 
support to the actual free volume space in the 
column was calculated and expressed as hold-up 
of the dispersed phase. 


By analogy with the results of an investigation 
by Srermour [7], Gayier et al. correlated the 
mean rate of rise v, of the droplets of the discon- 
tinuous phase relative to the moving continuous 
phase by the following expression : 


V V 
~=%(1—2) (1) 


e(1 —2) 


where 0, was defined as the limiting mean droplet 
velocity at zero continuous phase or the charac- 
teristic velocity. They indicated that this 
characteristic velocity could be determined by 


plotting 
(r+ 


against 2(1—2z) on arithmetic co-ordinates. 
Straight lines of slope ed, passing through the 
origin should be obtained. When data of the 
present investigation are plotted in this manner 
straightline correlationships result for specific 
packing and liquid systems giving the following 
basic correlation 


(2) 
where m — et». Fig. 2 shows such a plot for the 


system  kerosene-acetone-water using gin. 
Raschig rings as the packing. The data of the 


/| 


Kerosene acetone — woter 
nRaschig rings 


005 O10 O16 O20 O25 030 
X(I-X) 


Fic, 2. Typical plot of hold-up data. 


present investigation covering more than 300 
experimental runs were first analysed by plotting 


— vs. 2 (1 —2) 


and the values of m were determined by the 
method of least squares. These values of m are 
given in Table 3 along with the calculated values 
of m from the reported data of GayLer et al. [6] 
covering a wide range of tower diameters, liquid 
systems and different packings. 


400 600 ae 1000 1200 
ey 4] 
Fic. 3. Correlation of hold-up data. 
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Table 3. 


T. Srramamayya and G. S. 


Caleulated values for correlation of hold-up data 


System 


Rerosene 
acetone water 


Butyl acetate 
acetone water 


Benzene 
acetone water 


Methy! isobutyl 
ketone water* 


Buty! 
acetate water® 


Tower 


dia. Packing 


acetate water*® 


Iso-octane water*® 


Dibuty! 


ketone water* 


Ieo-octane 20°, 
glveerol® 


lin, R 
iin, B 


lin, KR 
tiin, 
lin, B 


jin, R 
fin, 
lin, RK 
jin, B 
lin, B 


jin, R 
jim. KR 
jim 
tin, 
jin, B 
lin, B 


jim, 
jin, B 


| jin, 
jin R 
tin, 
lin, B 
iin 


jin, 
tin, 
tin, B 


| 


a 


Ap 


WTS 


7200 


408-6 


Tt 
wee 


7174 


The dependence of m on physical properties of 
the liquid systems and packing characteristics 
should now be evaluated. 
of viscosity and interfacial tension it may be 
assumed that m should vary with g, @/e, p, and 


Raschig rings 
Berl sacidles 
Lessing rings 


» Calculated values from the data of al (6) 


Neglecting the effect 


12 
(= 


Ap. The following relationship results by the 
application of dimensional analysis ; 


(3) 


Tower | 
Benzene water® 
2 jin, 873-2 265 
jin. L 
jin. 
2 jin, BO 237 a 
2 jin v4 245-5 VOL 
2 jim, 75 13 
2 jin, 225-6 96 
| 
oe 
jim 468 
jin, B 52 uno 278 
tin. B 30s 
6 237 
lin, 7650 573 450-6 
| 
3 jm 4626 426 
jin, R 4 wiz 
lin, KR 7242 
K 
L 
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Hold-up in packed liquid liquid extraction columus 


The most useful value of n was found to be | from 
the analysis of experimental data, The validity 
of the assumptions made may now be examined 
by plotting the calculated values of vs. 
[(a/e* g) (p./Ap} on arithmetic co-ordinates. 
This is shown in Fig. 3 where the calculated 
values of m from the experimental data as well 
as from the data reported by Gayuen et al. [6] 
for a wide variety of liquid systems, packing 
types and sizes and a wide range of tower dia- 
meters are plotted against [(a eg) (p./ Ap) | 1/2 
It is observed that the data are satisfactorily 
correlated by the equation 


(4) 


m 0.683 | 


a p.) 
eg Ap | 
Average percentage deviation of the points from 
correlation is 10-07 per cent, Fifty-five per cent 
of the points lie within + 10 per cent of the corre- 
lation and 80 per cent of the points fall within 
the range of 15 per cent. 

The generalized correlations for hold-up in 
packed liquid-liquid extraction columns ranging 
from 2in, to 12im. in diameter packed with 
different sizes cf Raschig rings, Berl saddles and 
Lessing rings and covering a wide variety of 


liquid system may be obtained by substituting 
the value of m from equation (4) into equation 
(2). The following correlation results : 


a 

4 . 0-683 (1 — x) 
1 g Ap 

which may be regarded as useful in predicting 

hold-up below flooding in liquid-liquid extraction 

columns. 
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NOTATION 


superficial area of packing ft? ft? 
fractional voids in packed tower 

gravitational acceleration ft hr? 
rate of rise of droplets of discontinuous phase ft hr 
characteristic velocity ft hr 
superticial velocity of continuous phase (empty 
tower) ft hr 
superticial velocity of discontinuous phase (empty 
tower) ft hr 
fractional hold-up of dispersed phase in column for 
space 

density of continuous phase 
density of discontinuous phase 
density difference 


Ib 
Ib ft? 
Ib ft® 
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Book Reviews 


Gas Chromatography Abstracts 1958. Edited by 
BE. Knapman, Butterworths, London 1960, pp. 262. 
& @ 

Tur Gas Chromatography Discussion Group associated 
with the Institute of Petroleum was founded in 1955 and 
since then has been responsible for the organisation of 
regular bi-annual informal symposia and three biennial 
formal meetings of an international character. It is for 
these reasons that it is best known. However, among its 
other and less publicized activities it sponsors an abstract- 
ing service. The need for this has arisen through the quite 
remarkable acceleration in the rate of publication. At the 
end of 1959 something over 2100 papers had been published 
while at the present time between thirty and forty papers 
a monch are appearing, often in journals of limited avail- 
ability. These figures speak for themselves. 

The volume under review contains abstracts of 1468 
publications which comprise the literature from the 
effective birth of the subject in 1952 up to the end of 1958. 
Since they were never intended to serve as effective sub- 
stitutes for the original papers, they are, in most cases, 
very brief, and do little more than indicate the contents 
of the work discussed. Allied with the abstracts, however, 
is an extremely comprehensive index in which the reference 
numbers of the abstracts have been grouped under such 


headings as, for example, “ Applications and Specialized 
Separations,” “Sample Type” or “ Related Methods 
and Techniques.” Each of these general classifications is 
further broken down into more specific groupings. Thus, 
the location of those abstracts relevent to a particular 
problem is quick and simple. While the volume is of 
value to all it is evidently particularly useful to those not 
constantly active in the field and whose interest is in know- 
ing what can be done rather than how it is done, since 
by judicious use of the index and abstracts the amount 
of reading of original papers can be reduced to a minimum, 

There is no doubt that Dr. Knapman and his associates 
have eminently succeeded in their aim, and this is the 
more estimable since the time lapse from the inception 
of the abstracting service to the appearance of the volume 
was little over twelve months, In fact the group has been 
so effective that very shortly after the appearance of the 
volume under review that for 1959 became available, 
The standard of production of the volume is high and very 
few errors are detectable. The compilers and publishers 
are to be congratulated on the appearance of this most 
useful book, and it is to be hoped that future volumes in 
the series will maintain the standards set here. 


J. H. 
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Letters to the Editors 


Corrections to the paper, ‘‘ Perforated-plate performance "’ [1] 


Taste 3, shown below, gives corrected values for the 
weep points discussed in the reference paper. Table 3 
(p. 32), as published (Reference 1, p. 32) gives the co- 
ordinates of Point A for plate inserts 1 through 5, and 
gives the co-ordinates of Point B for plate inserts 6 and 
7. Weepage between Points A and B is generally not 
severe, but Point B is the generally accepted weep point. 
The new data points when compared to the line of 


HucumarK and O'Connett [2] give about as good 
agreement as the previously shown Figure 15 in the 
original article. Therefore, the remarks about the weepage 
correlation need not be changed in the original article. 


There is one other correction that should be noted in 
Figure 16, the Oscillation-point correlation, the last value 
of D,, should be 332 in. instead of 9/32 in. 


Table 3. Weep points 


3-68 
4-06 
346 
3-86 
+14 
3-46 
3-00 


2610 


E EE 


Department of Chemical Engineering, Lamar State College of Technology 
Beaumont, Texas 


R. A. McALuisTer 
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13 insert (Ib H,O/ (in. H,0) (Ib air/ F, 
96 number ft®) ft®) 
I 850 2-29 1275 148 
2-50 1220 
2000 2.69 1240 
350 2-33 1330 
1600 2.04 1260 19-6 
3000 2-65 1100 12-7 
i 350 2-18 1470 170 
1600 2-48 1280 148 
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IV 350 — 
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3000 2-53 1240 144 
v 350 2-50 1200 13-9 
1600 2-76 1050 121 
3000 2-95 1030 19 
vl 350 2-75 144 
1600 3-20 15-0 
3000 3-55 15-6 
vu 450 374 | 16-4 
720 16-1 
960 17-7 
1380 15-0 
1660 16-4 
2020 17-5 
2450 14-6 
— | 16-5 
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